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2  Physical  Parameter  Metrics 
for  Terrain 


A  terrain  physical  parameter  metric  is  a  set  of  constants  that  relates  the 
radiance  statistics  of  an  IR  image  to  the  physical  parameters  that  describe  the 
terrain.  Two  physical  parameter  metrics  are  required  for  each  terrain  feature 
because  both  the  standard  deviation  and  correlation  length  of  the  radiance  must 
be  related  to  the  physical  parameters  of  the  terrain  surface. 


Bare  Soil  Surface 

The  physical  parameters  necessary  to  describe  the  average  and  statistical 
properties  of  terrain  radiance  are; 

=  long  wave  emissivity  of  soil 

a^  =  short  wave  absorptivity  of  soil 

kj  =  heat  conductivity  of  soil 

Ps  =  thermal  diffusivity  of  soil 

Ys  =  soil  moisture 

m^  =  average  terrain  slope  of  soil  terrain 

These  parameters  can  be  written  as  a  physical  parameter  row  matrix  for  soil 
terrain  areas 

Ps  =  (1  Es  K  K  % 


where  the  unity  symbol  is  added  for  mathematical  convenience.  These  are  the 
minimum  number  of  physical  parameters  that  are  required  to  determine  terrain 
radiance. 
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Associated  with  the  row  matrix  of  physical  parameters  are  two  additional 
row  matrices  that  are  required  for  the  calculation  of  terrain  radiance  texture. 
These  are  the  physical  parameter  metric  matrices,  one  associated  with  the 
standard  deviation  of  terrain  radiance  and  the  other  associated  with  the  correla¬ 
tion  length  of  the  terrain  radiance.  The  physical  parameter  metric  associated 
with  the  standard  deviation  of  the  soil  terrain  radiance  is  as  follows: 


bs  =  (bo  bj  bj  b^  b^  b,) 


(2) 


where  the  b’s  are  empirical  constants  that  characterize  the  standard  deviation 
of  the  radiance  for  soil  terrain.  The  physical  parameter  metric  associated  with 
the  correlation  length  of  the  radiance  of  soil  terrain  is  as  follows: 


(^0  ^2  ^3  ^4  ^5  ^6^ 


(3) 


where  the  c’s  are  a  set  of  seven  empirical  constants  that  characterize  the 
correlation  length  of  the  radiance  of  the  bare  soil  surface. 


Grass  and  Short  Vegetation 

For  a  vegetated  surface  such  as  grass  or  secondary  vegetation  such  as 
bushes,  the  main  physical  parameters  that  are  necessary  for  the  radiance  cal¬ 
culation  are: 

Eiv  =  long  wave  emissivity  of  green  vegetation 

aiv  =  short  wave  absorptivity  for  green  vegetation 

mj^,  =  average  terrain  slope  for  green  vegetation 

=  long  wave  emissivity  of  brown  vegetation 

a^v  =  short  wave  absorptivity  for  brown  vegetation 

m^j^,  =  average  terrain  slope  for  brown  vegetation 

A  reasonable  approximation  might  be  m^  =  mi^  =  m^^,.  The  physical  para¬ 
meters  for  vegetation  can  be  written  as  two  row  matrices  of  the  form: 


Piv  =  (1  Eiv 


(4) 
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Corresponding  to  these  vegetation  parameters  are  the  physical  parameter  tex¬ 
ture  metrics  that  are  associated  with  the  standard  deviation  and  the  correlation 
length  of  the  radiance  texture  of  the  green  and  brown  vegetation  components. 
Thus  the  physical  parameter  metrics  associated  with  the  standard  deviation  of 
the  radiance  of  the  green  and  brown  vegetation  are: 


=  fro  *2  fr) 

fo.  -  (fo  f,  4  4)  (7) 


where  the  e’s  and  f’s  are  empirical  constants  that  describe  the  standard  devi¬ 
ation  of  the  vegetation  radiance.  The  physical  parameter  metrics  for  the  corre¬ 
lation  length  of  the  radiance  of  green  and  brown  vegetation  are: 


g, v  =  (go  gl  g2  gs)  (8) 

h, ,  =  (ho  h,  h,  ho)  (9) 


and  where  the  g’s  and  the  h’s  are  constants  that  describe  the  correlation  length 
of  the  radiance  signature  of  the  vegetation. 

Forest  Canopies 

For  forested  canopies,  the  basic  physical  parameters  that  are  necessary  for 
the  radiance  statistics  calculations  are: 

Ef  =  long  wave  emissivity  of  leaves 

Of  =  long  wave  absorptivity  of  leaves 

=  long  wave  emissivity  of  soil 

ttj  =  long  wave  absorptivity  of  soil 

Bf  =  short  wave  absoiptivity  of  leaves 
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If  =  leaf  area  index 

Ij  =  leaf  slope  distribution 

Tf  =  stomatal  resistance  to  water  vapor  diffusion 

where  heat  storage  has  been  ignored.  These  physical  parameters  for  forest 
canopies  are  represented  by  the  following  row  matrix: 


Pf  (1  £f  ctf  (x^  aj  If  Ij  tf) 


(10) 


The  corresponding  physical  parameter  texture  metric  that  describes  the  stan¬ 
dard  deviation  of  the  forest  canopy  radiance  is  described  by 


Pf  (Po  Pi  P2  P3  P4  P5  Pe  P7  Pg^ 


(11) 


while  the  texture  metric  that  describes  the  correlation  length  of  the  radiance  of 
forest  canopies  is  given  by; 


Qf  =  (Qo  Qi  Qs  04  Os  Oa  Ov  Og) 


(12) 


where  the  p’s  and  the  q’s  are  empirical  constants. 

The  physical  parameter  metrics  that  are  introduced  in  this  section  will  now 
be  used  to  develop  a  physics-based  radiance  texture  model. 
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3  Physics-Based  IR  Terrain 
Radiance  Texture  Modei 


The  physics-based  radiance  texture  model  is  loosely  based  on  the  numerical 
conclusions  of  some  previous  work  in  this  area  (Ben-Yosef  et  al.  1985,  1987; 
Ben  Yosef,  Wilner,  and  Abitbol  1986).  These  conclusions  are  that  (a)  the 
standard  deviation  of  the  terrain  radiance  is  proportional  to  the  solar  loading, 
and  (b)  the  correlation  length  is  inversely  proportional  to  the  solar  loading. 
These  conclusions  need  to  be  altered  as  can  be  seen  from  the  extreme  case 
where  the  emissivity  of  terrain  surface  is  taken  to  have  a  zero  value,  because 
in  this  case  the  thermal  emission  is  zero  irrespective  of  the  solar  loading.  A 
little  thought  suggests  that  the  more  appropriate  statements  should  be  (c)  the 
standard  deviation  of  the  terrain  radiance  is  related  to  the  average  terrain  radi¬ 
ance,  and  (d)  the  correlation  length  of  the  terrain  radiance  is  inversely  related 
to  the  average  terrain  radiance.  This  is  a  more  reasonable  relationship  because 
now  if  the  terrain  emissivity  has  a  zero  value,  the  mean  radiance  has  a  zero 
value,  the  standard  deviation  of  the  radiance  has  a  zero  value,  and  the  corre¬ 
lation  length  is  infinite. 


Terrain  Radiance  Texture  Modei  for  Bare  Soiis 

In  this  report,  the  standard  deviation  and  correlation  length  of  the  radiance 
for  bare  soil  areas  are  described  using  the  following  matrix  equations 


^Ns 


(13) 


p^,  =  rT^^cP 

r  Ns  s  s  s 


(14) 


where 
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Ons  =  Standard  deviation  of  soil  terrain  radiance 
r  =  resolution,  pixel/m 

Nj  =  average  radiance  of  bare  soil  in  a  specified  wave  band,  Wm‘‘  sf’ 

bs  =  empirical  physical  parameter  metric  for  the  standard  deviation  of  bare 
soil  radiance  as  represented  by  Equation  2 

P  s  =  transpose  column  physical  parameter  matrix  for  bare  soil  that  is 
obtained  from  Equation  1 

Pns  =  correlation  length  of  soil  radiance,  m 

Cs  =  empirical  physical  parameter  metric  for  the  correlation  length  of  soil 
radiance  as  represented  in  Equation  3 

6  =  empirical  exponent 

P  =  empirical  exponent 

r|  =  empirical  exponent 

K  =  empirical  exponent 

Equations  13  and  14  can  be  rewritten  as: 

On.  -  (15) 

P».  =  (15) 

where 

=  b„  +  b.e,  +  b-a  +  b,k  +  b.u,  +  b^^v  +  b,m^  (17) 

OS  0  Is  2s  3s  4“  s  5  *s  os 

Mps  =  Cq  +  CjE^  +  c^a^  +  Cjk  +  c^p^  +  c^y^  +  c^m^  (18) 

The  soil  terrain  radiance  is  given  by: 


10 
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r 

J 


N.  =  e.  Nj  dv 


where  NJ  =  N^Cv.TJ  =  Planck  function  for  bare  soil,  and  =  bare  soil  tem¬ 
perature  as  predicted  by  the  WESTHERM.  Thus  in  general: 

T  =  T  (e  .a^.k  ,p  ,Y  ,m  ,  weather,  time,  season)  (20) 


N*  =  N  te,,a.,k  ,p^,v,m^,Av,  weather,  time,  season) 


where  Av  =  Vj  -  Vj  =  wave  band.  Both  and  are  predicted  by  the 
WESTHERM.  Equations  15-18  show  that,  in  general,  16  constants  are 
required  for  the  semitheoretical  description  of  the  radiance  texture  for  a  soil 
terrain  area.  . 

A  simplification  of  Equations  15-19  can  be  obtained  by  considering  only 
zeroth  order  terms  which  gives: 


m;  =  b„ 


K  =  Co 


Then  Equations  15  and  16  become: 


=  b,r^N; 


Pns  = 


where  the  empirical  exponents  6,  p,  q,  and  k  are  given  in  Items  14  through  16 
of  Appendix  A.  Within  this  approximation,  only  two  constants,  b^,  and  c^,, 
need  to  be  determined  empirically  for  a  semiempirical  texture  model  for  soil 
terrain  radiance. 
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Terrain  Radiance  Texture  Model  for  Grass  and 
Short  Vegetation 


Consider  now  the  case  of  grass  or  short  vegetation  with  a  partial  cover  over 
soil.  The  standard  deviation  and  correlation  length  of  the  radiance  of  grass  or 
short  vegetation  covering  soil  is  given  by: 


®Nvs  ^v)  ®Ns  ^!v)  ^Ndv  ^iv  ®Nlv 


=  r^  {(1  -  X,) 


(1  -xj  I<.v  N.f]} 


Pnvs  =  In  {(1  -  x^)  exp  (-1/pj,,)  +  x^  [(1  -  x,^) 


exp  +  x,^  exp  (-1/Pn,v)]} 


where  o^s  ^nd  p^^  are- given  by  Equations  24  and  25,  and  where 


«^Ndv  =  f.v  Pdv  =  Nat 


On,  =  r<  e,  P,  =  r-^M^,  N* 


P»..  =  '■“Nd".  h„  P„  =  N,' 


Pn,.  =  r'N,;  g,.  P„  =  r-M  .  N,.' 
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where  the  empirical  exponents  8,  p,  ti,  and  k  are  given  in  Items  14  through  16 
of  Appendix  A 

where 

Cmvs  =  standard  deviation  of  radiance  of  vegetation  and  soil  terrain, 

Wm'^  sr' 

=  standard  deviation  of  radiance  for  brown  vegetation,  Wm'^  sr  ’ 

o^iv  =  standard  deviation  of  radiance  for  green  vegetation,  Wm'^  sr’' 

Pnvs  =  correlation  length  of  radiance  of  vegetation  and  soil  terrain,  m 

PNdv  =  correlation  length  of  radiance  for  brown  vegetation,  m 

Pfjiy  =  correlation  length  of  radiance  for  green  vegetation,  m 

x^  =  fraction  of  terrain  that  is  covered  by  vegetation 

Xjv  =  fraction  of  vegetation  cover  that  is  green 

Ndv  =  average  radiance  for  brown  vegetation  in  a  specified  wave  band, 
Wm’'^  sr' 

fdv  =  physical  terrain  metric  for  standard  deviation  of  radiance  from  dead 
vegetation  as  represented  by  Equation  7 

P*d„  =  transpose  column  physical  parameter  matrix  for  brown  vegetation 
that  is  obtained  from  Equation  5 

Njy  =  average  radiance  for  green  vegetation  in  a  specified  wave  band, 
Wm’^  sr' 

e,v  =  physical  terrain  metric  for  standard  deviation  of  radiance  from  green 
vegetation  as  represented  by  Equation  6 

=  transpose  column  physical  parameter  matnx  for  green  vegetation 
that  is  obtained  from  Equation  4 

hdv  =  physical  terrain  metric  for  the  correlation  length  of  the  radiance 
from  brown  vegetation  as  represented  in  Equation  9 

gj^  =  physical  terrain  metric  for  the  correlation  length  of  the  radiance  of 
green  vegetation  as  represented  by  Equation  8 

where  and  Mp^  are  given  by  Equations  17  and  18,  respectively,  and  where: 
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N„  -  e,,  J  N“  dv  (37) 

where 

N^'’  =  N^'’  (v,  T(j^,)  =  Planck  function  for  brown  vegetation 
Njj  =  Ni'’  (v,  Ti^,)  =  Planck  function  for  green  vegetation 
and  Tjv  and  T,^  =  temperature  of  brown  and  green  vegetation,  respectively. 
Therefore: 

Ndv  =  (Edv,  adv,  m^y,  Av,  weather,  time,  season) 

Niv  =  (e,„,  aiv,  m,y,  Av,  weather,  time,  season) 

Tdv  Tdv  (e^v,  a<jv,  m^^,  weather,  time,  season) 

Tiv  =  T],,  (e,^,  aiv,  mi^,  weather,  time,  season) 
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where  Av  =  Vj  -  v,  =  wave  band.  The  WESTHERM  computer  program 
calculates  T,y,  N^v.  and  N,.,.  Equations  26-35  show  that  16  constants  are 
required  to  predict  the  radiance  texture  of  grass  and  short  vegetation. 

If  no  bare  soil  spots  are  visible  and  the  vegetation  cover  is  complete,  then 
=  1  and  Equations  26  and  27  become: 


(1  OnjJv  +  X|^  C^;iv 


1/2 


(38) 


=  r'l 


(1  -  X,  )  Ml  +  X,, 


11/2 


Pnv  =  -In  f(l  -  x,^)  exp  (-l/Pi,J  +  x,^  exp  (-l/p^iv. 


(39) 


If  seasonal  effects  are  neglected  and  all  vegetation  is  considered  alive  in  the 
vegetation  cover,  then  Xj^  =  1  and  Equations  26  and  27  become  in  general 


o, 


Nlvs 


(1  -  xj  Ons  + 


_2 

^Nlv 


1/2 


(40) 


= 


(1  -  xj  m;  n“  +  X,  Ml  N, 


j2b 

ctIv  ^  Iv 


11/2 


Pnivs  =  -In  [(1  -  X^)  exp  (-l/p^,,)  +  x^  exp  (-l/p^iv)] 


(41) 


If  seasonal  effects  are  neglected  for  the  case  of  a  complete  vegetation  cover, 
then  Xy  =  1  and  x,^,  =  1  and  Equations  38  and  39  give: 


=  <Jniv  =  M^,, 


(42) 


Nv 


Pniv  =  rTM,:  Mp,, 


(43) 


A  simplification  occurs  in  Equations  26-43  by  neglecting  higher  order  terms 
and  considering  only  the  zeroth  order  terms  in  Equations  32-35  with  the  result 
that: 


Chapter  3  Physics-Based  IR  Terrain  Radiance  Texture  Model 


15 


ML  =  f 


m;,v  =  e„ 


ML  = 


Ml  =  go 


Within  this  approximation,  only  four  constants  are  required  to  describe  the 
radiance  statistics  of  vegetation.  Combining  Equations  22,  23,  26,  27  and 
44-47  gives: 


jL  =  {(1  -  xj  bo^  Nf  +  X,  [(1  -  X,,)  f„^ 


2 


2  XT  26 


and  pSvs  is  calculated  from  Equation  27  by  using  Equation  25  and 


pL  = 


pL  =  gof^i! 


with  the  result  that 


=  -In  {(1  -  xj  exp  (-1/pss) 


+  X  1(1  -  xj  exp  (-l/pNdv)+  X,  exp  (-I/Pnw) 
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which  are  valid  for  combined  vegetation  and  soil  terrain.  Equations  48-51 
show  that  within  the  zeroth  order  approximation,  only  six  constants,  b^,  c^,  e^, 
fg,  go,  and  ho,  are  required  to  be  determined  empirically  for  a  semitheoretical 
description  of  the  radiance  texture  for  a  vegetated  soil  terrain  surface. 


Terrain  Radiance  Texture  Modei  for  Forest 
Canopies 

The  standard  deviation  and  correlation  length  of  the  radiance  of  forest 
canopies  is  described  by: 


(52) 


pNf 


r’T^/  q,  P, 


(53) 


where  the  empirical  exponents  5,  p,  p,  and  k  for  forest  canopies  were  not 
determined  in  this  study,  and  where 

=  standard  deviation  of  forest  canopy  radiance 

Nf  =  average  radiance  of  forest  canopy  in  a  specified  wave  band,  Wm'^  sr  ' 

Pf  =  empirical  physical  parameter  metric  for  the  standard  deviation  of 
forest  canopy  radiance  as  represented  by  Equation  1 1 

=  transpose  column  physical  parameter  matrix  for  forest  canopies  that  is 
obtained  from  Equation  10 

=  correlation  length  of  forest  canopy  radiance,  m 

qj  =  empirical  physical  parameter  metric  for  the  correlation  length  of 
forest  canopy  radiance  as  represented  by  Equation  12 

Equations  52  and  53  can  be  rewritten  as: 
o  =  r^N/  M 

^Nf  f  of 
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(55) 


P»,  - 


where 


(56) 


=  Po  +  PlEf  +  P2«f  +  Ps^s  +  P4as  Ps^f  +  Pfilf  ^  PtIs 


(57) 


Mpf  =  Qo  +  q,et  +  q2af  +  qs^s 


+  q^a^  +  +  q,lj  +  q^l,  +  q^r^ 


The  average  forest  canopy  radiance  is  given  by; 

Nf  =  EjJ  nJ  dv  (58) 

v, 


where 

N^=  Nv  (v,  Tf)  =  Planck  function  for  forest  canopy 

and  where  Tf  =  forest  canopy  surface  temperature  as  predicted  by 
WESTHERM.  Therefore,  in  general: 


Tf  =  Tf(ef,af,lf,Gf,mf,  weather,  time)  (59) 


Nf  =  Nf(ef,af,lf,Gf,mf,Av,  weather,  time)  (60) 


where  Av  =  V2  -  Vj  =  wave  band.  Both  Tf  and  Nf  are  calculated  by  the 
WESTHERM.  Equations  52-57  show  that  for  the  general  case,  20  empirical 
constants  are  required  for  a  semiempirical  description  of  the  radiance  texture  of 
forest  canopies. 

Equations  52-57  can  be  simplified  by  dropping  all  first  order  terms  of  the 
forest  canopy  parameters  which  gives: 
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=  p, 


(61) 


M 


pf  =  qo 


(62) 


so  that  Equations  52  and  53  become: 


a 


0 

Nf 


=  r  W' 


(63) 


P^f  = 


(64) 


Therefore,  within  this  approximation  only  two  empirical  constants,  p^  and  q^, 
need  to  be  determined  to  determine  the  radiance  statistics  of  forest  canopies. 
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4  Terrain  Temperature 
Statistics 


Average  Radiance  and  Brightness  Temperature 

The  following  is  a  brief  summary  of  the  standard  theory  of  thermal  radia¬ 
tion  (Bramson  1968;  Hudson  1969).  The  radiance  in  a  wave  band  of  a 
surface  having  an  emissivity  value  e  and  a  physical  temperature  T  can  be 
obtained  from  the  following  approximate  form  of  Planck’s  law  (Weiss  and 
Scoggins  1989;  Weiss,  Scoggins,  and  Meeker  1992) 


N  =  eCf  M  III 


““ 

f 

exp 

C2 

- 1 

X  T  , 

- 

- 

(65) 


where  =  weighted  average  wavelength  of  the  thermal  radiation  in  the  wave 
band  A^,  and  where 


Cf  =  1.192  X  10*  W  pm^  m-2  sr'* 


(66) 


=  1.4338  X  1(T  pm  K 


(67) 


Equivalently,  the  radiance  can  be  written  in  terms  of  the  brightness  tempera¬ 
ture  Tg  as  follows 


N  =  Cf  AX  Xll 


exp 

f  C-^'  ^ 

-  1 

X  T„ 

- 

^  av  Bj 

- 

(68) 
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Using  Equations  65  and  68  gives  the  following  standard  results  (Bramson 
1968) 


e  = 


exp 


X  T 

^  av  ^ 


-  1 


exp 


X  T„ 

^  av 


-  1 


(69) 


^a'v 

Cl  AX  Xl^  N  -') 


p 

f  N  ^ 

-1 

In  ' 

1  +  e‘ 

exp 

Co  ’ 

-  1 

> 

X  T 

- 

(  -  J 

- 

(70) 


N  A  -1 


T  = 


C,  X 


In  (l 


e  cr  AX  X: 


N' 


Q"'  Kl 


'  1  +  e 

exp 

{ 

u 

_ 

-  1 

> 

- 

- 

^  av 

- 

(71) 


For  the  IR  wave  bands  of  interest  in  this  report,  namely  3  to  5  and  8  to  12  pm, 
the  following  Wien’s  law  approximations  to  Equations  65,  68,  and  69-71  are 
valid 
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( 


N  ~  e  Cf  AX,  Xjt  exp  ■ 

V 


X  T 


/ 

-  Cf'  AX  Xa^  exp  ■ 

V 


X 


av 


(72) 


In  e  ~  X3„^ 


1 

T 


.1 


(73) 


Tg  -  III  In-'  (cr  AX  X;'  N-') 


Q" 


In  e 


(74) 


-  T 


1  + 


Q" 


In  e 


T  ~  x;'  In-'  (e  cf  AX  x;^  N  -') 


~  T 

-n 


1  + 


T„  X. 


VI 


B  V 


In  e 


(75) 


~  T„ 


1  - 


T,  >... 


In  e 


Equations  72-75  are  only  approximations. 

As  an  example  of  the  use  of  Equations  65-75,  consider  the  8-  to  12-pm 
wave  band.  For  this  case,  AX  =  4  pm,  X^^  =  10.76  pm  and  therefore 
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cf  III  =  1332.5  K 


(76) 


Cr  AX  =  3305.8  Wm  sr 


-2  c.,--! 


(77) 


so  that  Equation  68  gives  the  radiance  for  the  8-  to  12-pm  wave  band  as 


N  -  3305.8  [exp  (1332.5/Tb)  -  1 


(78) 


3305.8  exp  (-  1332.5/Tb) 


Now  consider  the  example  of  the  special  case  of  the  3-  to  5-pm  wave  band  for 
which  AX,  =  2  pm  and  =  4.8  pm  so  that  for  this  case 


q""  X;„'  =  2987.1  K 


(79) 


cr  Xll  AX  =  93562.1  Wm'^  sr"' 


(80) 


and  therefore  the  radiance  for  the  3-  to  5-pm  wave  band  is  given  by 


N  =  93562.1  [exp  (2987.  l/T^) 


-  1 


(81) 


93562.1  exp  (-  2987.  ITT^) 


The  equations  given  above  are  useful  expressions  for  the  radiance  of  natural 
and  man-made  surfaces. 


Radiance  and  Temperature  Statistics 

This  section  determines  the  relationship  between  the  statistical  parameters 
of  the  terrain  radiance  and  the  statistical  parameters  of  brightness  temperature 
and  physical  temperature  of  the  terrain.  The  determination  of  the  standard 
deviation  of  the  brightness  temperature  in  terms  of  the  standard  deviation  of 
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the  radiance  is  accomplished  by  using  Planck’s  radiation  law,  which  is  written 
in  the  form  of  Equation  68.  The  terrain  radiance  calculation  that  includes  the 
effects  of  texture  can  be  obtained  from  Equation  68  as 


N  +  Of,  =  Cf  AX  X, 


-5  s 


exp 


X  (T„  +  Otb) 

av^  B 


-  1 


(82) 


where 

On  =  standard  deviation  of  the  terrain  radiance 

Oto  =  standard  deviation  of  the  terrain  brightness  temperature 

Then  a  Taylor  series  expansion  of  Equation  82  about  the  temperature  T^  keep¬ 
ing  only  the  linear  term  gives  the  standard  deviation  of  the  radiance  as 


O.J.B 


Q 


X  T„ 

V  V 


exp 


^  T„ 
V 


-  1 


-2 


^  -N  ^ 


Kl  Tb'^  N  o^b  exp 


Q 


^  A 


exp 


Q 


,^a.TB^ 


-  1 


(83) 


-  (cr  AX)-^  Xl  o.^,  (1  +  cr  AX  x;:  N-*)  in^ 


(1  +  cr  AX  x:^  N-‘) 


However,  in  the  wave  bands  3  to  5  and  8  to  12  pm  and  for  Tg  ~  300  K,  the 
following  Wien  approximation  is  valid 


exp 


X  T„ 
V 


))  1 


or  Cf  AX  Xll  N  -‘  ))  1 


(84) 


and  it  follows  approximately  that 
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~  Cf  Q  AX  Xay  Tg  exp  — — _  1 


cf  AX  c  Tb^  a^B  exp 


X  T„ 

av  B 


Q^'  Tb^  n  o^b 


N  Irf  (cr  A?1  ?.3t  N-')  o. 


Equations  83-87  can  be  used  to  calculate  the  standard  deviation  of  the  bright¬ 
ness  temperature  Ojb  terms  of  the  standard  deviation  of  the  radiance.  For 
example  Equation  83  gives 


=  (cr  A^)"'  Tb^  o^.  exp  - 


_  exp -  ”  1 

X  ^ 

av  B  j  I  av  B  i 


K  Tb  N->  o,.  exp  I  -  I  exp  |  -  1 


=  cf  Q"'  A^  C  N-"  o^.  (1  +  cf  AX  x:l  N-')-'  In-^ 


(1  +  C^'  A^  N-') 


while  the  approximate  Equation  85  gives 


Chapter  4  Terrain  Temperature  Statistics 


25 


-  (Cf  AXf  Xl 


~  (cr  Q"''  A^)’’  exp 

~  (c^y  X^^  T^  N-'  a, 

-  C^"-'  ?i:;  N  -'  ln-2  (cr  A?i  ^3^  N  y 

For  the  8-  to  12-pm  wave  band,  Equation  87  becomes 

~  2.27  X  10-'  Tb  [exp  (1332.5713)  -  l]  a^. 
~  2.27  X  10-'  Tb  exp  (1332.5713) 

~  7.50  X  10^  Tb  N-' 

~  1332.5  N-’  In-'  (3305.8/n)  o^. 

while  for  the  3-  to  5-pm  wave  band,  Equation  87  gives 
0,3  -  3.58  X  10-^  Tg  [exp  (2987. ITTg)  -  l]  <5^, 

~  3.58  X  10-’  Tb  exp  (2987.1713)  o^, 

~  3.35  X  10^  Tg  N  Oj, 

-  2987.1  N-‘  In-'  (93562.  1/n)  o^, 


(87) 


(88) 


(89) 
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These  equations  establish  the  relationship  between  the  standard  deviation  of  the 
radiance  and  the  standard  deviation  of  the  brightness  temperature. 

The  standard  deviation  of  the  radiance  can  also  be  expressed  in  terms  of  the 
standard  deviation  of  the  emissivity  and  the  standard  deviation  of  the  physical 
temperature  by  using  Equation  65  which  gives 


N  +  =  (e  +  o^)  Cf  III  jexp 


cT 

(T  +  O^) 


A  Taylor  series  expansion  of  Equation  90  keeping  only  the  linear  terms  gives 
the  following  results 


which  can  also  be  written  as 


=  Ne'^o^  +  Q'"  X,,'v  T"^  No.^  exp 


=  Ne-'o,  +  (e  Cf  A?.)-'  (1  +  e  Cf  A?i  K'  N'^) 


In^  (1  +  E  Cf  A?i  Xll  N"')  (-92^ 

In  a  manner  similar  to  that  shown  in  Equation  84,  the  following  inequalities 
are  valid  for  the  3-  to  5-  and  8-  to  12-pm  wave  bands 
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so  that  the  Wien  approximations  to  Equations  91  and  92  are 


Equations  91-94  relate  On,  Oj  and,  o,. 

Equations  83-85  and  91-94  can  be  used  to  estimate  the  standard  deviation 
of  the  emissivity.  A  comparison  of  Equations  83  and  92  gives 


=  (cr  Alf  K.  N  [a^  (l  +  cr  AX  N->)  In^  (l  +  Cf  AX  X,t  N ->) 
-  e-'o^  (l  +  e  AX  X^t  N  -')  Irf 

(l  +eCf  AXX3IN-')] 
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Combining  the  Wien  approximations  of  Equations  84  and  93  with  Equation  95 
gives 


oJ^  ~  C2  Kl  (Tb^  <^tb  ■  <^r) 


In^  (cr  AX  x:l  /V-') 


-  Oj.  In^  (e  AX, 


(96) 


Equation  96  can  be  used  to  estimate  but  requires  that  measurements  of 
be  available,  which  are  generally  difficult  to  obtain.  A  more  useful  expression 
can  be  obtained  from  Equation  96  by  assuming  the  approximation 

o,p  ~  o,rB  (97) 


Combining  Equations  96  and  97  gives 


oje  ~  Kl  -T-^)o,3 


(C,^')”'  X^^  [in^  (cr  AX  Kl  N-‘)  -  Itf  (e  Cf  AX  Kl  N-‘)J 


j\-l 


In^ 


e  +  2  In  e  In  (cf  AX  X.;^  N  *) 


-  2  In  e  (C^"']'’  In  (cf  AX  Xll  N'*)  0,^3 


(98) 


Combining  Equations  74  and  98  gives 


Og  —  2e  In  £  Tg'*  c,[.b 


(99) 
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Equation  99  shows  that  for  natural  terrain  materials  with  Tg  ~  300K,  o-rg  -  IK, 
and  a  range  of  values  of  emissivity  given  by  0. 1  <  e  <  0.98,  the  correspond¬ 
ing  range  of  is  given  by 

6.8  X  10-^  <  Oe  <  1-5  X  10'^ 


which  is  a  spread  of  only  a  factor  of  two. 

Finally,  the  relationship  between  the  correlation  lengths  of  the  radiance,  the 
physical  temperature,  and  the  brightness  temperature  is  given  by 

Pn  ~  Pt  ~  Ptb 

where 

=  correlation  length  of  the  terrain  radiance 
Pt  =  correlation  length  of  the  physical  temperature 
pjg  =  correlation  length  of  the  brightness  temperature 
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5  Determination  of  Empiricai 
Constants  of  the  Model 


Empirical  Description  of  Texture 


The  computer  model  WESTHERM  uses  heat  budget  algorithms  to  predict 
the  radiances  of  various  types  of  materials  usually  found  in  an  outdoor  thermal 
scene.  Types  of  materials  might  be  as  follows:  bare  soil,  green  vegetation, 
brown  vegetation,  bare  rock,  forest  canopy,  etc.  Each  of  these  material  types 
can  have  many  varieties.  For  example,  bare  soil  may  have  many  different 
compositions,  with  different  thermal  properties.  Similarly,  green  vegetation 
may  consist  of  grass,  tall  weeds,  or  bushes,  each  with  different  thermal  proper¬ 
ties.  Forest  canopies  may  be  of  decidous  trees  or  evergreen  trees,  each  with 
different  thermal  properties,  depending  upon  the  season  of  the  year.  The  pur¬ 
pose  of  algorithms,  such  as  WESTHERM,  is  to  predict  the  radiance  of  a  given 
material  over  a  diurnal  cycle,  given  the  geographical  location  and  the  day  of 
the  year.  The  problem  is  that  these  programs  predict  only  the  average 
radiances  of  the  material  under  a  given  set  of  circumstances.  In  the  real  world, 
there  is  a  natural  fluctuation  about  the  average — the  radiance  of  a  material  is 
said  to  have  texture.  This  texture  is  customarily  measured  in  terms  of  two 
texture  parameters — standard  deviation  and  correlation  length. 

The  standard  deviation  is  the  usual  statistical  parameter  computed  about  the 
mean  value  of  the  radiance  over  a  region  composed  of  the  same  material  type. 
The  correlation  length  is  defined  in  terms  of  the  autocovariance  function  (the 
Fourier  transform  of  the  power  spectrum)  of  the  radiance  over  a  region  com¬ 
posed  of  the  same  material  type.  This  is  assumed  to  be  identical  to  the  linear 
autocorrelation  coefficient  (also  called  Pearson’s  r).  The  linear  autocorrelation 
coefficient  of  a  series  Xj  is  defined  by  the  equation 


Y,  {x-x)(x^_-x) 


(102) 


The  function  f(t)  is  assumed  to  be  exponential  of  the  form  y  =  exp  (-t/p), 
where  p  is  the  correlation  length  of  the  series  Xj.  Strictly  speaking,  this  is  not 
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true,  since  F  has  a  range  from  -1  to  1,  and  the  exponential  function  is  always 
positive.  So  it  is  assumed  that  F(t)  can  be  approximated  by  the  exponential 
function  on  the  interval  0  <  t  <  p. 


Power  Law  Model 

The  physics-based  model  of  the  IR  radiance  statistics  for  terrain  includes 
several  constants  that  have  to  be  determined  empirically  from  measured  IR 
images  of  terrain  areas  that  include  bare  soil,  grass,  and  forested  canopies  for 
which  the  physical  parameters  and  local  weather  conditions  are  known.  A 
statistical  analysis  of  these  IR  images  yields  the  average,  standard  deviation, 
and  correlation  length  of  the  radiance  corresponding  to  known  time  of  day, 
terrain  characteristics,  and  weather  conditions.  The  physics-based  radiance 
statistics  model  is  then  fit  to  these  statistical  parameters,  and  the  unknown 
constants  of  the  model  are  determined. 

The  assumption  upon  which  the  texture  parameter  model  is  based  is  that 
both  the  standard  deviation  and  the  correlation  length  of  radiance  can  be 
approximated  by  a  power  law.  That  is,  if  N  represents  radiance,  and  r  is  the 
camera  resolution,  and  o^,  and  pj^:  represent  standard  deviation  and  correlation 
length  of  radiance,  respectively,  then  it  is  assumed  that  these  quantities  can  be 
approximated  by  the  equations 

(103) 


and 


p^  =  (104) 

where  b,c  5  and  (3  are  model  constants  which  must  be  determined  by  a  statisti¬ 
cal  analysis  of  measured  data,  and  k  and  r]  are  empirically  determined  con¬ 
stants.  The  current  working  hypothesis  is  that  k  =  -1  and  fi  =  1.  This  model 
uses  the  normalized  radiance  obtained  by  dividing  the  radiance  by  the  average 
radiance  over  a  diurnal  cycle.  The  normalized  radiance  is  given  by 

N 

n  =  —  (105) 


where  is  the  average  radiance  over  a  diurnal  cycle.  In  the  same  manner, 
the  normalized  resolution  is  written  as 
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r 


(106) 


where  r|.gf  a  reference  resolution  which  in  this  report  is  chosen  to  have  the 
value  rj.gf  =1.0  pixel/m.  In  this  way,  the  coefficients  b,  5,  c  and  p  refer  to 
normalized  radiances  and  normalized  resolutions  so  that  the  dimensions  of 
b  and  c  are  given  by 

[b]  = 'Wm~'^sr~^  [c]  =  m 

A  collection  of  computer  programs,  referred  to  collectively  as  WESTEX 
have  been  developed  to  empirically  determine  the  model  constants  and  then 
predict  the  diurnal  variation  of  the  terrain  texture  through  the  use  of 
WESTHERM.  A  listing  of  the  WESTEX  modules  appears  in  Appendix  A. 

Yuma,  AZ,  Data 

In  the  spring  of  1993,  data  were  collected  at  Yuma,  AZ,  in  the  form  of  IR 
images  of  desert  scenes.  These  images  were  taken  once  or  twice  an  hour  over 
a  period  of  days.  Figures  2-7  give  the  measured  diurnal  variations  of  the 
average  radiance,  standard  deviation  of  the  radiance,  and  correlation  length  of 
the  radiance.  These  data  are  obtained  from  measured  IR  images  of  terrain  that 
were  acquired -at  the  Yuma,  AZ,  test  site  in  the  8-  to  l2-pm  wave  band  for 
soil  areas  with  soil  and  grass  mixtures  (Rivera  1994a).  Two  of  these  diurnal 
cycles  were  selected  for  analysis.  The  diurnal  cycles  covered  the  time  period 
from  midtright  to  midnight  for  the  dates  April  8  and  April  26,  1993.  The 
terrain  at  the  Yuma  test  area  consists  of  several  identifiable  materials,  includ¬ 
ing  desert  pavement  and  mixed  grass/bare  soil.  For  each  material,  and  for 
each  image,  the  average  radiance,  standard  deviation  of  radiance,  and  correla¬ 
tion  length  of  radiance  were  computed.  Figures  8-13  show  the  average  radi¬ 
ance,  standard  deviation  of  the  radiance,  and  correlation  length  of  the  radiance 
as  a  function  of  time  for  the  2  days  of  measured  data  that  were  used  in  this 
report.  Figures  14  and  15  show  the  measured  values  of  the  standard  deviation 
of  the  radiance  and  the  correlation  length  of  the  radiance  versus  the  measured 
values  of  the  average  radiance  for  2  days  of  data  collection,  while  Figures  16 
and  17  show  the  measured  values  of  the  standard  deviation  and  correlation 
length  versus  the  normalized  values  of  the  average  radiance.  Figures  14-17 
show  that  the  power  law  assumptions  given  by  Equations  103  and  104  are  not 
valid  for  the  fuU  range  of  the  average  radiance  that  occurs  during  a  diurnal 
cycle. 

It  was  found  necessary  to  divide  the  day  into  seven  time  zones,  and  the 
measured  standard  deviations  and  correlation  lengths  of  the  radiance  were  fit  in 
each  zone  by  the  method  of  least  squares  to  the  measured  normalized  radiances 
using  Equations  103  and  104  and  using  the  data  represented  in  Figures  16  and 
17.  The  results  appear  in  Figures  18-21,  which  show  the  least  square  fits  of 
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LWB.  26APR93:  DESERT  PAVEMENT  MEAN  TEMPERATURE  BOUNDED  BY  STD.  DEV.  AND  CORR.  LENGTH 


00:00  02:00  04:00  06:00  08:00  10:00  12:00  14:00  16:00  18:00  20:00  22:00  24:00 


Measured  values  of  the  diurnal  variation  of  the  average,  standard 
deviation,  and  correlation  length  of  the  radiance  for  desert  pave¬ 
ment  at  Yuma  Test  Site  on  24  March  1993 


LWS.  26APR93:  GRASSABARE  SOIL  MEAN  TEMPERATURE  BOUNDED  BY  STD.  DEV.  AND  CORR.  LENGTH 


00.00  02:0C  04:00  06:00  08:00  10:00  12:00  14:00  16:00  18:00  20:00  22:C0  24:00 


Measured  values  of  the  diurnal  variation  of  the  average,  standard 
deviation,  and  correlation  length  of  the  radiance  for  grass  and  bare 
soil  at  Yuma  Test  Site  on  24  March  1993 
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LWB,  24MAR93:  DESERT  PAVEMENT  MEAN  TEMPERATURE  BOUNDED  BY  STD.  DEV.  AND  CORR.  LENGTH 


Measured  values  of  the  diurnal  variation  of  the  average,  standard 
deviation,  and  correlation  length  of  the  radiance  for  desert  pave¬ 
ment  at  Yuma  Test  Site  on  8  April  1993 


Figure  5.  Measured  values  of  the  diurnal  variation  of  the  average,  standard 
deviation,  and  correlation  length  of  the  radiance  for  grass  and  bare 
soil  at  Yuma  Test  Site  on  8  April  1993 
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Measured  values  of  the  diurnal  variation  of  the  average,  standard 
deviation,  and  correlation  length  of  the  radiance  for  desert  pave¬ 
ment  at  Yuma  Test  Site  on  26  April  1993 


Measured  values  of  the  diurnal  variation  of  the  average,  standard 
deviation,  and  correlation  length  of  the  radiance  for  grass  and  bare 
soil  at  Yuma  Test  Site  on  26  April  1993 


Chapter  5  Determination  of  Empirical  Constants  of  the  Model 


Chapter  5  Determination  of  Empirical  Constants  of  the  Model 


37 


Chapter  5  Determination  of  Empirical  Constants  of  the  Model 


Chapter  5  Determination  of  Empirical  Constants  of  the  Model 


□□ 


□ 


Chapter  5  Determination  of  Empirical  Constants  of  the  Model 


Chapter  5  Determination  of  Empirical  Constants  of  the  Model 


41 


Figure  12.  Measured  diurnal  variation  of  the  standard  deviation  of  the  radiance  for  desert  pavement  at  Yuma,  AZ,  on  26  April  1993 
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not  adequate  to  describe  the  24-hr  diurnal  variation 
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A  display  of  the  measured  correlation  length  of  the  radiance  versus  the  measured  radiance  for  desert  pavement  at  Yuma,  AZ, 
on  8  April  and  26  April  1993.  The  solid  and  dotted  lines  represent  a  power  law  fit  to  the  data  and  show  that  the  24-hr  diurnal 
variation  is  not  describable  by  this  means 
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Figure  16.  A  comparison  of  the  measured  standard  deviation  of  the  radiance  and  the  normalized  measured  radiance  tor  desert  pavement 
at  Yuma,  AZ,  on  8  April  and  26  April  1993.  The  solid  and  dotted  lines  are  best  tit  power  law  tits  to  the  data  and  are  presented 
to  show  the  inadequacy  of  this  description  of  the  tuil  24-hr  diurnal  variation 
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Figure  18.  A  2-day  least  squares  parameter  fit  using  a  power  law  for  the  standard  deviation  of  the  radiance  versus  the  radiance  in  seven 
radiance  bands  is  compared  with  the  measured  data  from  8  April  1993  at  Yuma,  AZ 
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Figure  19.  The  2-day  least  squares  parameter  fit  using  a  power  law  of  the  correlation  length  of  the  radiance  versus  the  radiance  over 
seven  radiance  bands  is  compared  with  the  measured  data  for  8  April  1993  at  Yuma,  AZ 
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Figure  20.  A  2-day  least  squares  parameter  fit  using  a  power  law  to  represent  the  standard  deviation  of  the  radiance  versus  the  radiance 
in  each  of  the  seven  radiance  bands  is  compared  with  the  measured  data  for  26  April  1993  at  Yuma,  AZ 


Figure  21 .  The  2-day  least  squares  parameter  fit  using  a  power  iaw  to  represent  the  correlation  length  of  the  radiance  versus  the  radiance 
in  seven  radiance  bands  is  compared  with  the  measured  data  for  26  April  1993  at  Yuma,  AZ 


the  standard  deviation  and  correlation  length  of  the  radiance  versus  the  average 
radiance  for  the  power  laws  given  by  Equations  103  and  104.  In  this  way,  the 
coefficients  b,  6,  c,  and  P  are  determined.  The  power  law  representations  of 
the  measured  standard  deviation  and  correlation  length  of  the  radiance  are 
displayed  versus  time  in  Figures  22-25  in  which  the  solid  lines  refer  to  the 
least  square  fit  that  appears  in  Figures  18-21.  Each  parameter  is  plotted  first 
versus  radiance,  since  that  is  the  domain  in  which  the  power  function  fit  is 
made,  and  secondly  the  parameter  is  plotted  versus  time  of  day,  so  that  the 
diurnal  variation  may  be  seen  in  the  time  domain.  These  curves  were  gener¬ 
ated  using  the  actual  measured  values  of  radiance,  rather  than  the  values  of 
radiance  predicted  by  WESTHERM.  The  model  constants  were  determined  by 
taking  the  weighted  averages  of  the  model  constants  for  the  2  days  according 
to  the  equations 


b  =  exp 


In  6,  +  n2  In 


c  =  exp 


^  1  .  ^ 

rtj  In  Cj  +  rt2  In 


n.  +  rto 

V  *  ^  y 


6  = 


«,6j  +  rt252 


n^  +  n^ 


(108) 


(109) 


(110) 


P  = 


/ll  +  «2 


(111) 


where  n,  and  n2  represent  samples  sizes  for  the  2  days.  Only  the  model  con¬ 
stants  for  desert  pavement  were  used  in  the  development  and  verification  of 
the  model. 


Data  Reduction  Procedure 


Given  a  set  of  measured  data  {(Xi,yi),  1  <  i  <  n},  it  is  desired  to  find  con¬ 
stants  a  and  b  such  that  if  y*  =  axj^,  then  y*  will  give  the  best  approximation 
to  y  in  the  least  squares  sense.  In  actual  practice,  however,  rather  than  fit  the 
power  curve  directly,  what  is  actually  fit  is  the  natural  logarithm  of  the  power 
cu^e:  In  (y^*)  =  In  (a)-i-  b  In  (Xj).  IJ  =  In  (Xj),  Yj  =  In  (y^,  and  Y*  -  In 
(yj*),  then  this  equation  becomes  Y,*  =  A  -i-  BXj,  where  A  =  In  (a)  and  B  =  b. 
The  values  of  A  and  B  can  be  found  as  the  constants  of  the  ordinaiy'  regres¬ 
sion  line  for  the  measured  values  of  (Xi,Yj).  Then  a  =  exp  (A)  and  b  =  B. 

The  diurnal  cycle  was  divided  into  seven  segments,  and  model  constants 
were  determined  for  each  of  the  seven  segments  for  standard  deviation  and 
correlation  length.  This  was  accomplished  using  program  Tit’  (a  compilation 
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Figure  23.  Comparison  of  the  least  squares  fit  of  the  correlation  length  of  the  radiance  versus  the  radiance  with  the  measured  correiation 
length  of  the  radiance  expressed  as  a  function  of  time  for  8  April  1993 


of  the  radiance  versus  rad 
993  at  YUMA,  AZ 


FigurG  25.  A  comparison  of  the  least  squares  fit  of  the  correlation  length  of  the  radiance  versus  the  radiance  with  the  measured  values  of 
the  correlation  length  of  the  radiance  when  expressed  as  function  for  time  for  26  April  1993  at  Yuma,  AZ 


of  files  fitc.c  and  fitf.f  in  Appendix  A),  which  used  data  files  ss040893.dat  and 
ss042693.dat  in  Appendix  B  (Once  Tit’  has  been  compiled,  one  types  ’fit 
ss040893.dat’  to  run  the  program).  The  output  of  this  program  consists  of  data 
files  of  the  type  ss040893.prms  and  ss042693.prms  in  Appendix  B.  These 
files  contain  the  model  constants  b,c,5,  and  (3,  as  well  as  information  needed 
when  averaging  the  two  data  sets.  The  values  of  k  and  r|  are  needed  for  mak¬ 
ing  adjustments  based  upon  camera  resolution.  The  program  ’avg’  (formed  by 
compiling  files  avgc.c  and  avgf  f  listed  in  Appendix  A)  averages  the  data  from 
several  output  files  of  the  type  ss(M0893.prms  and  ss042693.prms,  and 
produces  a  file  of  the  sort  ss.avg  in  Appendix  B.  To  produce  an  input  file  for 
avg  one  must  concatenate,  into  a  single  file,  all  the  files  with  extension 
’prms’  which  one  wishes  to  average,  taking  care  to  remove  the  extra  lines 
defining  k  and  Tj,  leaving  only  the  first  two  such  lines.  This  produces  a  file  of 
the  type  ss.all  in  Appendix  B.  One  runs  the  compiled  program  by  typing  ’avg 
ss.aU.’  This  produces  the  file  ’ss.avg.’ 

The  full  model  assumes  that  model  constants  b,c,6,  and  P  for  each  of  the 
seven  segments  (called  bands),  are  known  for  four  different  materials:  bare 
soil,  green  vegetation,  brown  vegetation,  and  forest  canopy.  To  have  a  full 
working  model,  one  needs  to  have  a  model  database  of  the  sort  represented  by 
the  file  modelYDP.db  in  Appendix  B.  This  is  done  by  concatenating  files 
such  as  ss.avg  into  a  single  file  and  removing  the  extra  lines  defining  k  and  Tj. 
Thus,  one  would  need  real  world  diurnal  data  of  the  type  found  in  files 
ss040893.dat  and  ss042693.dat  for  each  of  these  four  material  types.  No  such 
data  are  presently  available.  For  bare  soil,  data  from  desert  pavement  near 
Yuma,  AZ,  was  used.  This  material  is  probably  at  the  extreme  end  of  the 
range  of  what  would  be  called  bare  soil.  Since  there  are  no  data  for  the  other 
three  categories,  model  YDP.db  contains  zeros  for  the  model  constants  in  those 
categories.  Once  one  has  a  model. db  file,  this  serves  as  the  input  file  for  the 
program  ’block’  (created  by  compiling  files  blockc.c  and  blockf  f).  Program 
’block’  turns  the  file  model.db  into  a  Fortran  block  data  subroutine,  called 
block.f  This  routine  serves  as  one  of  the  modules  for  the  texture  prediction 
program  ’texgen’  (compiled  from  files  texgen.c  sdaclor.f  fsdaclor.f  and  block.f 
found  in  Appendix  A). 


Selection  of  Diurnal  Bands 


Since  the  model  predicts  the  values  of  standard  deviation  and  correlation 
length  in  terms  of  the  value  of  radiance,  some  method  was  needed  to  divide 
the  diurnal  cycle  (which  is  a  time  cycle)  according  to  the  value  of  the  radiance 
since  radiance,  and  not  time,  is  the  independent  parameter  of  the  model. 
Dividing  the  day  into  bands  determined  by  time  of  day  produced  unsatisfactory 
results.  During  a  diurnal  cycle,  the  radiance  cycles  through  a  range  of  values, 
going  from  a  minimum  value  (typically  around  sunrise)  to  a  maximum  value 
during  the  heat  of  the  day,  then  back  to  a  lower  value  after  sunset.  The  seven 
bands  were  determined  according  to  the  value  of  the  minimum,  average,  and 
maximum  radiance  during  a  diurnal  cycle.  First,  all  radiances  in  the  measured 
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data  cycle  are  normalized  by  dividing  by  the  average  radiance  over  the  cycle. 
Let  us  denote  the  minimum  and  maximum  normalized  radiances  by  and 
n^ax>  respectively.  In  addition  to  these  values,  arbitrary  values  s  and  t  are 
selected  and  adjusted  empirically.  Experimentation  showed  that  values  of  s  = 
0.5  and  t  =  0.84  produced  good  results.  In  general,  s  and  t  must  satisfy  0  <  s 
<  1  and  0  <  t  <  1.  The  values  of  s  and  t  are  used  to  determine  two  more 
boundaries  on  the  normalized  radiance  axis,  nj  and  nj,  with  n^  =  (1  -  sjn^jn  +  s 
and  Hj  =  (1  - 1)  +  tn^ax-  s  is  the  proportionate  distance  of  n^  between 
rijnin  and  1,  while  t  is  the  proportionate  distance  of  itj  between  1  and  n^^a^. 

Over  any  diurnal  cycle,  it  is  always  true  that  n^^jn  <  n^  <  1  <  i\  <  n^^ax’  Next, 
the  day  is  divided  into  two  time  zones:  Zone  1  is  that  period  of  the  day  before 
the  maximum  radiance  is  attained,  and  zone  2  is  the  period  after  the  maximum 
radiance  is  attained.  A  given  data  point  falls  into  one  of  the  seven  bands 
according  to  the  following  rules  based  upon  the  measured  value  n  of  the  nor¬ 
malized  radiance: 

Band  1:  zone  1,  and  n  <  Oj 
Band  2:  zone  1,  and  nj  <  n  <  1 
Band  3:  zone  1,  and  1  <  n  <  n^ 

Band  4:  zone  1,  and  ni  <  n  <  n^^^x 
Band  5:  zone  2,  and  1  <  n  <  n^j^x 
Band  6:  zone  2,  and  nj  <  n  <  1 
Band  7;  zone  2,  and  n  <  n^ 

These  bands  are  shown  in  Figures  18-21.  Notice  that  bands  1  and  7  are 
actually  merged  into  a  single  band.  Data  points  lying  in  bands  1  and  7  are 
merged  together  by  the  ’fit’  routine  to  find  common  model  constants  for  both 
bands.  Bands  2,  4,  and  6  are  transition  bands  whose  constants  are  computed 
in  such  a  way  as  to  avoid  functional  discontinuities  in  the  model.  So,  in  prac¬ 
tice,  only  data  points  in  bands  1,  3,  5,  and  7  are  actually  needed  to  determine 
the  model  constants  for  all  seven  bands.  Model  constants  for  bands  2,  4,  and 
6  are  computed  from  the  model  constants  of  bands  1,  3,  5,  and  7  by  the  ’fit’ 
routine.  However,  the  model  constants  for  these  transition  bands  are  recom¬ 
puted  at  run-time  since  they  depend  upon  the  actual  radiances  predicted  by 
WESTHERM.  The  equations  used  to  compute  the  model  constants  in  bands  2, 
4,  and  6  are  presented  below. 

Computation  of  Model  Constants  for  Bands  2,  4, 
and  6 

For  the  models  at  hand,  the  task  is  to  define  two  functions  that  are  piece- 
wise  power  functions  on  an  interval  [n^ij,>^ax]’  where 

<  1  <  n^  <  The  two  power  functions,  Pi(n)  and  PjCn),  are 
defined  as  follows: 
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‘  for  ^ 

for  n<n<\ 
a.^n  ^  for  \  <  n  <  n^ 
ay ‘for  n,<n< 


ajfi  ’  for  <  n  < 

P2^ri)=  '  ay^  for  <  n  <  1 

ay^for  I  <n<  ^ 

where,  for  some  s  and  t  satisfying  0  <  s  <  1  and  0  <  t  <  1, 


n,  =  (l-t)l+m 


By  the  nature  of  the  application,  the  values  of  aj,  aj,  aj,  ay,  bj,  53,  bj,  and  by 
are  given.  It  is  also  given  that  aj  =  ay,  and  bj  =  by.  The  problem  is  to  deter¬ 
mine  the  values  of  a2,  a^,  b2,  b4,  and  bg  which  will  make  Pj(n)  and  P2(n) 

continuous  on  the  interval  [r^i„,r^,J,  and  which  wiU  make  P,(n^,,)  =  P2(i\n,,). 
In  order  for  these  conditions  to  be  satisfied,  the  following  equations  must  hold: 


ay^  =  ay^ 


=  a^n, 
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(119) 


'^4^max  ~  '^S^max 


(120) 


h  b^  (121) 

From  Equations  117  and  120  it  follows  that  ^  =  a3>  2nd  a^  =  a^,  independen¬ 
tly  of  the  values  of  s,  t,  n^j^,  and  n^^^.  From  Equations  1 16  and  1 17,  we 
conclude  that 


^2  =  ^1 


In  -  In  03 
In 


From  Equations  118  and  119  we  conclude  that 


In 


64=. 


a-iti 


3'U 


"5"max 
7 - ^ 


In 


^”max  ^ 


and 


04=05/1, 


b5-b4 

max 


From  Equations  121  we  conclude  that 
In  -  In  ^5 


66=^1+- 


In  n. 


(122) 


(123) 


(124) 


(125) 


Model  Validation 

After  deriving  model  constants  for  standard  deviation  and  correlation  length 
for  the  seven  bands,  based  upon  the  Yuma  desert  pavement  data,  it  was  possi¬ 
ble  to  attempt  a  verification  of  the  model.  In  the  spring  of  1994,  data  were 
collected  at  Grayling,  MI  (Rivera  1994b).  This  was  not  diurnal  data,  so  a  full 
validation  could  not  be  conducted.  IR  images  from  several  days  were  ana¬ 
lyzed  in  order  to  determine  radiance,  standard  deviation  of  radiance,  and  corre¬ 
lation  length  of  radiance  from  a  patch  of  sandy  road.  While  there  was  no 
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reason  to  expect  that  the  model  parameters  for  desert  pavement  at  Yuma,  AZ, 
would  be  identical  to  model  parameters  for  a  sandy  test  track  at  Grayling,  MI, 
it  was  felt  that  they  should  be  sufficiently  similar  and  that  the  predicted  diurnal 
variation  would  produce  curves  of  the  expected  shape,  perhaps  being  off  by 
some  constant  multiple.  This  is,  in  fact,  what  was  found.  Only  three  or  four 
data  points  were  available  for  each  day  for  which  data  were  available.  Fig¬ 
ures  26-33  show  a  comparison  of  the  predicted  and  measured  radiances  and 
brightness  temperatures  for  3  days  of  measurements  at  the  Grayling  test  site. 
Actual  model  constants  for  the  Grayling  site  are  unknown  because  diurnal 
radiance  measurements  were  not  acquired.  In  modifying  the  Yuma  model 
constants,  only  the  values  of  b  and  c  were  changed.  The  values  of  the  expo¬ 
nents  p  and  6  were  unchanged.  The  altered  model  constants  are  found  in  the 
file  ’modelGTT.db’  in  Appendix  B.  Brightness  temperature,  Tg,  and  standard 
deviation  of  brightness  temperature,  Ojg,  were  computed  from  radiance,  N, 
and  standard  deviation  of  radiance,  Ojyf,  using  the  following  approximations  on 
the  8-  to  12-pm  range: 

j,  ^  1332.5 


OjB  -  2.27  X  10-Vj  [exp  (127) 

which  are  obtained  from  Equations  78  and  88.  Results  were  converted  to 
Celcius  using  the  relation 

/ft  =  Tfi  -  273.15  (128) 


^tb  ~  ^TB  (129) 

where  tj,  and  are  the  average  and  standard  deviations  of  the  brightness 
temperature  measured  in  degrees  Celsius. 


Use  of  WESTEX  Program 

In  its  present  form,  the  model  can  be  used  to  predict  diurnal  variation  in 
standard  deviation  of  radiance  and  correlation  length  of  radiance  for  either 
desert  pavement  at  Yuma,  AZ,  or  sandy  test  track  at  Grayling,  MI,  for  any 
date  for  which  a  met  file  exists.  The  met  file  is  a  prerequisite  for  running 
WESTHERM,  from  which  the  diurnal  radiance  is  predicted.  The  TEXGEN 
program  is  produced  by  compiling  the  files  ’texgen.c’,  ’sdaclor.f’,  ’fsdaclor.f , 
together  with  one  of  the  bdata.f  files,  either  ’bdataYDP.f  or  ’bdataGTT.f . 
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The  program  is  designed  to  be  run  in  a  Unix  environment  with  command  line 
parameters.  One  may  type  ’texgen’  to  get  a  list  of  the  required  parameters. 
They  are,  in  order:  (a)  xv,  the  percent  of  vegetation  versus  bare  soil  in  the 
surface  material,  (b)  xlv,  the  percent  of  the  vegetation  that  is  green, 

(c)  resolution,  the  camera  resolution  in  pixels  per  meter,  and  (d)  infile,  the 
name  of  the  output  file  from  WESTHERM  containing  the  diurnal  radiance 
profile  (such  as  the  file  08994.bare  in  Appendix  B);  outfile,  the  name  of  the 
file  you  wish  the  output  of  texgen  to  go  to. 

Since  there  are  presently  no  model  constants  available  for  green  vegetation, 
brown  vegetation,  or  forest  canopy,  xv  and  xlv  should  be  set  to  0.0.  Actually, 
the  model  can  be  run  in  this  mode  on  any  material  by  substituting  the  model 
parameters  for  Yuma  desert  pavement  with  the  model  parameters  for  the  mate¬ 
rial  in  question.  In  order  to  determine  those  model  parameters,  it  is  necessary 
to  have  diurnal  information  about  the  material  in  the  format  found  in,  for 
example,  the  file  ss040893.dat  in  Appendix  A.  If  the  information  is  available 
and  cannot  conveniently  be  put  into  that  format,  then  the  ’fetch’  subroutine  of 
the  ’fitf.f’  module  will  have  to  be  rewritten  to  accommodate  the  different  for¬ 
mat.  (This  is  a  trivial  exercise  for  a  programmer.)  Running  ’fit’  on  the  data 
wiU  produce  a  file  that  can  be  manipulated  according  to  the  procedure  outlined 
in  part  5,  Data  Reduction  Procedure,  to  produce  an  appropriate  bdata.f  file.  If 
this  procedure  is  used  to  find  model  parameters  for  green  vegetation,  brown 
vegetation,  and  forest  canopy,  then  the  model  parameters  may  be  cut  and 
pasted  into  the  model.db  file  to  produce  a  complete  bdata.f  module. 
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Figure  26.  Theoretical  prediction  of  the  time  dependence  of  the  radiance,  standard  deviation  of  the  radiance,  and  the  correlation  length  of 
the  radiance  compared  with  the  corresponding  measured  values  of  these  quantities  at  Grayling,  Ml,  on  30  March  1994 


Chapter  5  Determination  of  Empirical  Constants  of  the  Model 


63 


Figure  27.  Theoretical  prediction  of  the  time  dependence  of  the  brightness  temperature,  standard  deviation  of  the  brightness  temperature, 
and  the  correlation  length  of  the  brightness  temperature  compared  with  the  corresponding  measured  values  of  these  quantities 
at  Grayling,  Ml,  on  30  March  1994 
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Figure  31 .  Theoretical  prediction  of  the  time  variation  of  the  brightness  temperature,  standard  deviation  of  the  brightness  temperature,  and 
the  correlation  length  of  the  brightness  temperature  compared  with  the  measured  values  of  these  quantities  for  Grayling,  Ml,  on 
3  April  1994 
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Figure  33.  Theoretical  prediction  of  the  time  variation  of  the  brightness  temperature,  standard  deviation  of  the  brightness  temperature,  and 
the  correlation  length  for  the  brightness  temperature  compared  with  the  measured  values  of  these  quantities  for  Grayling,  Ml, 


6  Conclusions  and 
Recommendations 


Conclusions 


This  report  describes  a  physics-based  terrain  radiance  model  that  predicts 
the  standard  deviation  and  correlation  length  of  the  IR  radiance  of  terrain  areas 
such  as  soil,  grass  and  shrubs,  and  forest  canopies.  The  predicted  radiance 
statistics  are  expressed  in  terms  of  sensor  characteristics  such  as  resolution  and 
wave  band,  and  in  terms  of  environmental  conditions  such  as  weather,  time  of 
day,  and  the  physical  characteristics  of  terrain.  This  simple  model  requires  the 
determination  of  several  empirical  constants  by  fitting  the  model  to  measured 
IR  images  at  several  terrain  sites  including  soU,  grass,  and  forest  canopy  areas. 
AVhen  these  constants  are  determined,  the  model  wUl  have  the  capability  of 
predicting  the  radiance  statistics  for  any  time  of  day  and  for  any  terrain  area 
where  the  weather  conditions  and  terrain  characteristics  are  specified.  The 
terrain  texture  model  that  is  described  in  this  report  can  be  used  to  generate 
texture  for  uniform  areas  of  soil,  grass,  soil-grass  combination,  and  tree  cano¬ 
pies  that  often  appear  in  synthetically  generated  IR  scenesAmages. 

The  study  of  IR  terrain  texture  modeling  produced  the  following 
conclusions: 

a.  A  model  for  predicting  IR  terrain  texture  can  be  developed  by  using 
simple  power  law  equations  to  connect  the  average  values  of  the  stan¬ 
dard  deviation  and  correlation  length  of  the  radiance  in  a  uniform  area 
of  an  IR  image  to  the  average  values  of  the  point  radiance  in  the  image. 
Because  the  average  of  the  point  radiance  values  is  predicted  by 
WESTHERM  in  terms  of  terrain,  weather,  season,  and  time  of  day,  it 
follows  that  the  terrain  texture  model  can  predict  the  standard  deviation 
and  correlation  length  of  the  radiance  in  terms  of  these  environmental 
conditions. 

b.  A  set  of  site-independent  and  weather-independent  empirical  coeffi¬ 
cients  can  be  determined  and  used  to  relate  the  standard  deviation  and 
correlation  length  of  the  radiance  to  the  average  value  of  the  radiance 
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in  an  IR  image.  These  coefficents  are  determined  from  measured  IR 
images  of  uniform  terrain  areas. 

c.  Planck’s  radiation  law  can  be  used  to  relate  the  standard  deviation  and 
correlation  length  of  the  radiance  to  the  standard  deviation  and  correla¬ 
tion  length  of  the  brightness  temperature,  physical  temperature,  and 
terrain  surface  emissivity. 


Recommendations 

The  study  of  the  modeling  of  IR  texture  of  terrain  areas  suggests  the  fol¬ 
lowing  recommendations: 

a.  Diurnal  IR  image  data  should  be  acquired  and  processed  for  various 
terrain  types  in  addition  to  the  desert  pavement  considered  in  this 
report.  These  terrain  types  include  the  following:  bare  soils;  asphalt, 

.  gravel,  and  concrete  roads;  grass  and  tree  canopies.  These  data  should 
then  be  included  in  WESTEX  texture  prediction  capability. 

b.  The  terrain  texture  model  that  is  described  in  this  report  should  be  used 
in  conjunction  with  the  ARMA  texture  generation  algorithm  to  generate 
synthetic  IR  scenes/images  that  include  a  texture  component.  The 
synthetically  generated  scenes/images  should  be  compared  with  mea¬ 
sured  images  to  check  for  realism. 

c.  A  theoretical  study  should  be  initiated  to  examine  the  physical  basis  of 
IR  texture  and  to  use  measured  IR  data  to  determine  terrain  parameter 
variability. 
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Appendix  A 

Data  Reduction  Modules  and 
WESTEX  Texture  Prediction 
Modules 


Appendix  A  Data  ReductionAA/ESTEX  Texture  Prediction  Modules 


Item  A.l:  Contents  of  file  'fitc.c' 

#include  <stdio.h> 

#include  <s  tr ing . h> 
main  (argc,argv) 
int  argc ; 
char  *argv[ ] ; 

( 

if  (  argc<2  )  { 

puts("usage:  fit  filename"); 
exit(O) ; 

) 

fit_(argv[l] ) ; 
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Item  A. 2:  Contents  of  file  'fitf.f' 


subroutine  fit (infile) 


************************************************************************ 
* 

* 

* 

★ 

* 

* 

* 

* 

* 

* 


* 
* 
* 
* 
* 
* 
* 
* 
* 
* 

★*************★★*******************★******★***★*********'*•******•*•*★****** 


This  program  takes,  as  input,  a  file  containing  radiance  and 
standard  deviation  of  radiance  (in  watts  per  square  meter  per 
steradian)  and  correlation  length  (in  meters)  of  one  of  the 
four  material  types:  bare  soil,  green  short  vegetation,  brown 
short  vegetation  and  forest  canopy.  The  data  must  cover  a 
complete  twenty  four  hour  diurnal  cycle  with  at  least  one  data 
point  per  hour.  The  file  must  be  arranged  and  formatted 
to  be  read  by  the  fetch  routine. 


real  x(60) ,y(60) ,z(60) , times (60) ,xnorm(60) 

character*13  date (60) 

character*20  feature 

character*20  infile 

real  b(7) ,delta(7) ,c(7) ,beta(7) 

integer  num(4) 

real  eta, kappa 

data  eta  /l.O/ 

data  kappa  /-l.O/ 

* 

*  Read  in  the  data  file.  x  -*  radiance,  y  —  standard  deviation  of 

*  radiance,  and  z  =  correlation  length. 

* 

call  fetch(resolution,n,x,y ,z ,date, feature , infile) 

* 

*  Normalize  the  data  by  dividing  by  the  average  radiance  over 

*  the  diurnal  cycle,  xnorm  -  normalized  radiance. 

* 

call  no rmal(n,x, xnorm) 

* 

*  For  each  of  the  seven  segments,  do  a  least  squares  curve  fit  for 

*  standard  deviation  and  correlation  length  in  terms  of  normalized 

*  radiance  according  to  the  model  y  —  p^x'^q ,  where  p  and  q  are 

*  constants  whose  value  is  determined  by  the  method  of  least  squares 

* 

*  For  the  first  call,  find  b  and  delta  such  that  y  —  b*xnorm'^ delta 

*  for  each  of  the  seven  segments. 

* 

call  lsq(n,xnorm,y,b , delta, num) 

* 

*  For  the  second  call,  find  c  and  beta  such  that  z  -  c*xnorm'^beta 

*  for  each  of  the  seven  segments . 

* 

call  lsq(n,xnorm,z ,c ,beta,num) 

* 

*  Divide  the  b  and  c  coefficients  by  the  resolution  factors  to  unbias 
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*  the  coefficients.  In  the  model,  these  factors  will  be  remultiplied 

*  ^  resolution  factors  appropriate  to  the  cameras  being  simulated. 

*  The  resolution  factor  for  standard  deviation  is  resolution'^eta  and 

*  for  correlation  length,  it  is  resolution'kappa.  The  provisional 

*  values  are  eta  -  1  and  kappa  -  -1.  If  later  experimental  evidence 
points  to  different  values  for  eta  and  kappa,  appropriate  changes 

^  must  be  made  in  the  "data"  statements  at  the  top  of  this  subroutine. 

call  resolve ( resolution , eta , kappa , b , c) 
xavg  -«  xnorm(n+l) 

write(*, ' (a,f6.2) ' )  "Average  diurnal  radiance  -  ",xavg 
xmax  -«  xnorm(n+2) 
xmin  -  xnorm(n+3) 
rn  -  xmin/xavg 
rx  -  xmax/xavg 

call  results (date , feature , rn, rx , num, b . delta , c ,beta , infile 
&  , eta, kappa) 

end  !  fit 

subroutine  fetch(resolution . n . x , y , z , date , feature , filename) 

real  x(60) ,y(60) , 2 (60) 

character*13  date (60) 

character*72  heading(4) 

character*12  filename 

character*20  feature 

open  (UNIT-1 , STATUS-' old' , FILE-filename , ERR-40) 

read  (1, ' (a, 2f9 . 5) ' )  feature , rl , r2 
resolution  -  2.0/(rl  +  r2) 

do  loop-=l,4 

read  (l,'(a)')  heading(loop) 
end  do 


^read  d . ' (2 (a . lx) , lx , f8 . 3 , 5x , 2 ( f 6 . 4 , 7x) ) ' , END-30)  date (n) , feature , 
n-n+1  x(n),y(n).z(n)  ‘ 

go  to  20 
continue 
close(l) 
n  —  n  -  1 
return 


write  (*,'(a)')  "***  FILE  NOT  FOUND  ***" 
stop 

end  !  fetch 
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************************************************************************ 
subroutine  normal (n,x,xnorm) 
real  x(60) ,xnorm(60) 
character*12  infile 

sumx  -  0,0 
xmax  “  x(l) 
xmin  -  x(l) 
do  i-l,n 

sumx  -  sumx  +  x(i) 
if  (  xmax  .It.  x(i)  )  then 
xmax  -  x(i) 
end  if 

if  (  xmin  .gt.  x(i)  )  then 
xmin  -  x(i) 
end  if 
end  do 

xbar  -  sumx/real(n) 
xnorm(n+l)  -  xbar 
xnorm(n+2)  -=  xmax 
xnorm(n+3)  -=  xmin 
do  i«l,n 

xnorm(i)  -  x(i)/xbar 
end  do 
return 
end 

★****^****************************************************************** 
subroutine  resolve ( resolution , eta , kappa , b  , c ) 
real  b(7),c(7) 
real  eta, kappa 

do  i-=l ,  7 

b(i)  -  b(i)/resolution**eta 
c(i)  -  c(i)/resolution**kappa 
end  do 

return 

end  !  resolve 

subroutine  lsq(n,x,y , a,b ,num) 

real  x(60) ,y(60) 
real  a(7),b(7) 
integer  num(4) 
real  lx,ly 
data  s/0.5/ 
data  t/0.84/ 

xbar  -  x(n+l) 
xmax  -  x(n+2) 
xmin  -  x(n+3) 
rx  -  xmax/xbar 
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rn  -  xmin/xbar 
rb  -  1.0  -  t  +  t*rx 
ra  “  (1.0  -  s)*rn  +  s 
Asumlxly  -  0.0 
Asumlx  “  0.0 
Asumly  -=  0.0 
Asumlxlx  -  0.0 
Bsumlxly  -  0.0 
Bsumlx  -  0.0 
Bsumly  -0.0 
Bsumlxlx  -  0.0 
Csumlxly  -  0.0 
Csumlx  -0.0 
Csumly  -  0.0 
Csumlxlx  -  0.0 
Dsumlxly  -  0.0 
Dsumlx  -  0.0 
Dsumly  -  0.0 
Dsumlxlx  -  0.0 

do  i-1 , 4 

nuin(i)  -  0 
end  do 

kflag  -  0 

do  i  -  l,n 

lx  -  alog(x(i)) 

ly  -  alog(y(i)) 

if  (  x(i)  .It.  ra  )  then 

Asumlxly  -  Asumlxly  +  lx*ly 
Asumlx  -  Asumlx  +  lx 

Asumly  •=  Asumly  +  ly 

Asumlxlx  -  Asumlxlx  +  lx*lx 
num(l)  -  num(l)  +  1 

else  if  (  x(i)  . ge .  1.0  .and.  x(i)  .It.  rb 
.and.  kflag  .eq,  0  )  then 
Bsumlxly  -  Bsumlxly  +  lx*ly 
Bsumlx  -  Bsumlx  +  lx 

Bsumly  -  Bsumly  +  ly 

Bsumlxlx  -  Bsumlxlx  +  lx*lx 
num(2)  -  num(2)  +  1 

else  if  (  x(i)  . ge .  rb  . and.  kflag  . eq .  0  )  then 
Csumlxly  -  Csumlxly  +  lx*ly 
Csumlx  -  Csumlx  +  lx 

Csumly  -  Csumly  +  ly 

Csumlxlx  -  Csumlxlx  +  lx*lx 
num(3)  -  num(3)  +  1 
if  (  rx  -  x(i)  .le.  0.0001)  then 
kflag  -  1 
end  if 

else  if  (  x(i)  . ge .  1.0  .and.  kflag  . eq .  1  )  then 
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Dsumlxly  —  Dsumlxly  +  lx*ly 
Dsumlx  -  Dsumlx  +  lx 
Dsumly  —  Dsumly  +  ly 
Dsumlxlx  “  Dsumlxlx  +  lx*lx 
Tiuni(4)  -  num(4)  +  1 
end  if 

end  do 

p  -  real(nuin(l) ) 

b(l)  -  (p*Asumlxly  -  Asumlx*Asuinly)/(p*Asumlxlx  -  Asuinlx*Asxiinlx) 
a(l)  —  exp(  (Asumlxlx*Asumly-Asumlx*Asumlxly) 

&  /  (p*Asumlxlx-Asumlx*Asuunlx)  ) 

p  -  real(num(2)) 

b(3)  -  (p*Bsumlxly  -  Bsumlx*Bsumly)/(p*Bsumlxlx  -  Bsuiiilx*Bsumlx) 
a(3)  -  exp(  (Bsxiiiilxlx*Bsuinly-Bsumlx*Bsumlxly) 

&  /  (p*Bsumlxlx-Bsumlx*Bsumlx)  ) 

p  —  real(nuin(3) ) 

b(4)  -=  (p^Csumlxly  -  Csumlx*Csuinly)/(p*Csuinlxlx  -  Csumlx*Csuinlx) 
a(4)  “  exp(  (Csumlxlx*Csuinly-Csiiinlx*Cs\iinlxly) 

&  /  (p*Csumlxlx-Csuinlx*Csuinlx)  ) 

p  =  real(nuin(4) ) 

b(5)  •=  (p^Dsumlxly  -  Dsuinlx*Dsuinly)/(p*Dsuinlxlx  -  Dsunilx*Dsumlx) 
a (5)  -=  exp(  (Dsuinlxlx*Dsuinly-Dsumlx*Dsuinlxly) 

&  /  (p*Dsuinlxlx-Dsuinlx*Dsumlx)  ) 

★*** 

*  Make  the  correct  adjustments  for  a(2) , a(6) ,b(2)  and  b(6) 

*  to  render  the  functions  continuous. 

'krkrk'k 

a(2)  -=  a(3) 

b(2)  =  b(l)  +  (  alog(a(l))  -  alog(a(3))  )/alog(ra) 

a(7)  ^  a(l) 

b(7)  «  b(l) 

a(6)  =  a(5) 

b(6)  -  b(7)  +  (  alog(a(7))  -  alog(a(5))  )/alog(ra) 

return 

end  !  Isq 

************************************************************************* 
subroutine  results(date,  feature, rn,  rx,num,b, delta,  c , beta,  infile 
&  , eta, kappa) 

real  x(60) ,y(60) ,2(60) ,b(7) , delta(7) ,c(7) ,beta(7) 
integer  num(4) 
character*13  infile, out 
character*72  heading(2) 
real  eta, kappa 
character*13  date (60) 
character*7  day 
character*20  feature 
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*  Name  the  output  file  by  striping  the  extensions  off  the  input  file 

*  and  adding  the  extension  ".prms" 

* 

locdot  •*  index ( inf ile ) 

out  -  infile(l:locdot-l)  //  ".prms" 

heading(l)  -  "band  b  delta  c  beta  " 

heading(2)  -  " .  .  . n 

day  -  date(l) (1:7) 

open(unit-l , file-out) 

write(l, ' (a,f6.2) ')  "eta  -",eta 
write(l/(a,f6,2)')  "kappa  ,kappa 
write(l , ' (a, 2x, a) ' )  day, feature 

write(l,'(a,4(i3,2x))')  "Bandsize  1, 3,4.5:  "  .  (nuin(i)  ,  i-1 , 4) 
write(l,  (a,f5.3)')  "Minimum  Radiance  /  Average  Radiance  —  ",rn 
write(l , ' (a , f5 . 3) ' )  "Maximum  Radiance  /  Average  Radiance  -  ”,rx 
write(l, ' (a)')  heading(l) 
writed, '(a)')  heading(2) 
do  i-1, 7 

write  (l/(i3,2x,g,2x,f7.3,2x,f8.2,2x,f7,3)')  i, 

&  b(i) ,delta(i) ,c(i) ,beta(i) 

end  do 
close  (1) 

write(*,'(2a)')  "Output  file :  ",out 
return 

end  !  results 
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Item  A. 3:  Contents  of  file  'avgc.c 


#include  <stdio.h> 

#include  <string.h> 
main  (argc,argv) 
int  argc ; 
char  *arg;v[  ]  ; 

{ 

if  (  argc<2  )  { 

puts("usage:  avg  filename"); 
exit(O) ; 

} 

average_(argv[l] ) ; 

) 
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Item  A. 4:  Contents  of  file  'avgf.f 


subroutine  average (infile) 

* 

★ 

*  This  subroutine  takes,  as  input,  the  assembled  values  of  b,  delta,  * 

*  c  and  beta  from  several  different  days ,  and  averages  the  values .  * 

*********************************************^^^^^^^^^^^^^^^^^^^^^* 

real  b(10 , 7) , delta(10 , 7) , c(10 , 7) ,beta(10,7) 
real  m(10),rx(10) 
integer  sainpsize(10 ,4) 
character*20  infile 

call  f etcb(n, b , delta , c , beta , infile , samps ize , rn , rx) 
call  ave (n,b , delta , c ,beta , samps ize , rn, rx) 

call  store(n, b, delta, c, beta) 

end  !  main 

subroutine  fetch(n,b , delta, c , beta , inf ile , sampsize , rn  rx) 

real  b(10 , 7) , delta(10 , 7) , c(10 , 7) ,beta(10 , 7) 

real  eta, kappa 

integer  sampsize(10,4) 

real  rn(lO) , rx(lO) 

character*7  date 

character*8  variable 

character*18  string 

character*20  material 

character*20  infile , outf ile 

character*38  words 

character*48  heading(2) 

open  (unit-1,  file-inf ile, status-' old' ,err-100) 

loc  -  index (inf ile, ". ") 
if  (loc  .ne.  0)  then 

outfile  -  infile(l:loc-l)  //  ".avg" 
else 

write(*,'(a)')  "***INPUT  FILE  MUST  CONTAIN  EXTENSION***- 
STOP 
end  if 

open  (unit-2, file-outf ile) 

read  (1 , ' (a, f6 . 2) ' )  variable, eta 
write(2,'(a,f6,2)')  variable, eta 
read  (1 , ' (a, f6 , 2) ' )  variable , kappa 
write(2, ' (a,f6.2) ')  variable , kappa 
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20 


n  -  1 

read  (1 (2a) end-50)  date , material 
read  (1 , ' (a,4(i3 , 2x) ) ' )  string, (sampsize(n, i) , i-l ,4) 
read  (1, ' (a, f5 . 3) ' )  words, rn(n) 
read  (1 , ' (a, £5 . 3) ' )  words, rx(n) 
do  i-l, 2 

read(l , ' (a) ' )  heading(i) 
end  do 

do  i-l, 7 

read  (1, ' (i3,2x,g,2x,f7.3,2x,f8.2,2x,f7.3) ' )  j,b(n,i), 

&  delta(n, i) ,c(n, i) ,beta(n, i) 

write  (2,'(i3,2x,fll.8,2x,f7.3,2x,f8.2,2x,f7.3)')  j,b(n,i), 
6e  delta(n,  i)  ,c(n,  i)  ,beta(n,  i) 

end  do 
n  -  n+1 
go  to  20 

50  n  —  n-1 

close(l) 

write (2 , ' (2a) ' )  "AVERAGES :  " , material 

do  i-l, 2 

write(2 , ' (a) ' )  heading(i) 
end  do 

write (*,' (2a) ' )  "Output  file:  ",outfile 
return 

100  write  (*,'(3a)^)  "File  infile (1 : length) '  not  found," 
stop 

end  !  fetch 

***********★************************************************:*•***★★***** 
subroutine  ave(n,b , delta, c ,beta, samps ize , rn,rx) 

* 

*  This  subroutine  averages  the  values  of  b, delta, c  and  beta  for  each 

*  of  the  seven  segments  over  several  data  sets. 

* 

real  b(10 , 7) , delta(10 , 7) , c(10 , 7) ,beta(10 , 7) 

real  rn(lO) , rx(lO) 

integer  sampsize(10 ,4) 

integer  k(4) 

data  s  /O . 5/ 

data  t  /0.84/ 

k(l)  -  1 
k(2)  -  3 
k(3)  -  4 
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k(4)  -  5 


do  j“l,4 

suml  -  0.0 
sum2  “0.0 
sum3  “  0.0 
sum4  “  0.0 
size  “  0.0 

do  i-l,n 

sizenow  -  real ( samps ize ( i, j )) 
size  “  size  +  sizenow 

suml  -  suml  +  alog(  b(i.k(j))  )*sizenow 
sum2  -  sum2  +  delta(i,k(j))*slzenow 
sum3  -  sum3  +  alog(  c(i.k(J))  )*sizenow 
sum4  -  sum4  +  beta(l,k(j ) )*sizenow 
end  do 

b(n+l,k(j))  -  exp(  suml/size) 

delta(n+l ,k(j ) )  -  sum2/size 
c(n+l,k(j))  —  ®xp(  sum3/size  ) 

beta(n+l ,k(j ) )  -  sum4/size 

end  do 

b(n+l,7)  “  b(n+l,l) 
delta(n+l,7)  -  delta(n+l,l) 
c(n+l,7)  “  c(n+l,l) 
beta(n+l,7)  “  beta(n+l,l) 

rnsum  “0.0 
rxsum  “  0.0 

do  i“l , n 

rnsum  “  rnsum  +  rn(i) 
rxsum  “  rxsum  +  rx(i) 
end  do 

rmin  “  rnsum/real(n) 
rmax  -  rxsum/real(n) 
ra  “  (1.0  -  s)*rmin  +  s 
rb  “  1.0  -  t  +  t*rmax 

call  adj us t ( ra , rb , rmax , n , b , delta) 
call  adjust(ra,rb, rmax, n.c, beta) 

return 
end  !  ave 


subroutine  adj  us  t ( r a , rb , rmax , n , a , b ) 
real  a(10 , 7) ,b(10 , 7) 
c  “  1.0/alog(ra) 
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* 

* 


a(n+l,2)  -  a(n+l,3) 
a(n+l,6)  «  a(n+l,5) 

b(n+l,2)  -  b(n+l,l)  +  (  alog(a(n+l , 1) )  -  alog(a(n+l , 3) )  )/alog(ra) 
b(n+l,6)  -  b(n+l,7)  +  (  alog(a(iH-l , 7) )  -  alog(a(n+l , 5) )  )/alog(ra) 
b(n+l,4)  -  b(n+l,3)  +  (  alog(a(n+l , 3) )  -  alog(a(n+l , 4) )  )/alog(rb) 
a(n+l,4)  -=  a(n+l ,  5)*rmax**(  b(n+l,5)  -  b(n+l,4)  ) 


return 


end  !  adjust 


subroutine  store (n,b , delta, c , beta) 

real  b(10,7) , delta(10 , 7) , c(10 , 7) ,beta(10,7) 


Unit  2  is  opened  in  the  'fetch'  routine 


do  i«l,7 

write(2, '  (i3,2x,g,2x,f7.3,2x,f8.2,2x,f7.3) ' )  i, 
>  i)  I  delta(n+l ,  i)  ,  c(n+l ,  i)  ,beta(n+l ,  i) 

end  do 


close(2) 


return 
end  !  store 
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Item  A. 5:  Contents  of  file  'blockc.c' 


#include  <stdio.h> 

#include  <string.h> 
main  ( argc , argv) 
int  argc ; 
char  *argv[ ] ; 

{ 

if  (  argc<2  )  { 

puts("usage:  block  filename’*); 
exit(O) ; 

) 

block_(argv[l] ) ; 
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Item  A. 6:  Contents  of  file  'blockf.f' 

subroutine  block( inf ile) 

***********************************************************************^ 


*  * 

*  This  subroutine  takes  as  input  a  file  containing  the  values  of  eta  * 

*  and  kappa  and  the  values  of  b,  delta,  c  and  beta  for  all  five  * 

*  segments  and  for  each  of  the  material  types:  bare  soil,  green  vege-  * 

*  tation,  brown  vegetation,  and  forest  canopy.  These  values  are  then  * 

*  reorganized  in  the  form  of  a  FPRTRAN  block  data  module.  This  module  * 

*  is  compiled  with  the  texgen.c,  sdaclor.f  and  fsdaclor.f  modules  to  * 

*  build  the  texgen  module.  The  block  data  module  is  named  ”bdata.f".  * 

*  * 


*******************************************^*************************** 

real  b(7 ,4) , delta(7 ,4) . c(7 ,4) ,beta(7,4) 
real  eta, kappa 
character*20  infile 

call  fetch(b, delta, c, beta, eta, kappa, infile) 
call  makebd(b , delta , c , beta , eta , kappa) 

end 

************************************************************************ 
subroutine  fetch(b , delta, c , beta, eta, kappa, infile) 
real  b(7 ,4) , delta(7 ,4) , c(7 ,4) ,beta(7,4) 
real  eta, kappa 
character*40  line 
character*20  stuff 
character*7  variable 
character*20  infile 

open(unit=l ,  f ile*=inf  ile  ,  status*' old'  ,  err-=100) 
read(l , ' (a, f 6 . 2) ' )  variable, eta 
read(l , ' (a, f6 . 2) ' )  variable , kappa 
do  i*l,4 

read(l , ' (a) ' )  stuff 
read(l , ' (a) ' )  line 
read(l , ' (a) ' )  line 
do  j-1,7 

read(l. ' (i3 , 2x, g, 2x, f7 . 3 ,2x, f8 . 2 , 2x, f7 . 3) ' ) 

&  k,b( j , i) , delta( j , i) , c( j , i) ,beta(j , i) 

end  do 
end  do 
close(l) 
return 

100  write (*,' (2a) ' )  infile,"  not  found." 

stop 

end 
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subroutine  niakebd(b ,  delta ,  c , beta ,  eta , kappa) 
real  b(7 , 4) , delta(7 ,4) , c(7 ,4) ,beta(7,4) 
real  eta, kappa 

open(unit«l , f ile-"bdata . f " ) 

write(l, ' (6x,a)')  "block  data" 
write(l, ' (6x,a) ' ) 

&"conimon  /values/  b(7 , 9 ,4)  . c(7 . 9 .4) . delta(7 ,4)  .beta(7 .4)  . eta.kappa" 
write(l,'(6x,a)')"real  eta. kappa" 
write(l,'(6x,a)')"data  b" 

writed, '(5x,a,6(f9.5. a) )')"&/". b(l.l) ,",",b(2.1).*.",b(3,l),"." 

^  „  .l>(4,l).".",b(5,l),",-,b(6,’l)i-!- 

writed, ' (5x,a,f9.5,a,a)')"&".b(7,l),".","56*0.0." 
write (1, ' (5x,a.6(f9.5,a)) ' .b(l. 2) , " , " ,b(2 , 2) ! " , " ,b(3 , 2) , "  ,  " 

*  .b(4,2),”,",b(5.2),",",b(6,2),"  " 
write(l, '(5x.a,f9.5,a,a)')"&".b(7.2).","."56*0.0." 
write(l,'(5x,a,6(f9.5,a))')"&".b(1.3),".".b(2.3),’".".b(3.3),"." 

^  .b(4,3),",",b(5.3),".".b(6,3)  "  " 

writed,  '  (5x,a,f9.5,a.a) '  ,b(7.3)  ,","."56*0  0  " 

writed,  '(5x,a,6(f9.5,a))')"&",bd,4),",",b(2.4)i",".b(3.4),"," 

^  ,b(4,4).".",b(5,4).",".b(6  4)  "  " 

wrLte(l,'(5x,a.f9.5,a,a)')"&",b(7,4),",","56*0.0/" 

write(l, '(6x,a)')"data  c" 

writed.  '(5x,  a,  6(f9.5,  a)) ')"&/".  cd,l).",".c(2.  l),",",c(3.1),"," 

^  .c(4,l) ,","  .c(5,l) ’c(6,’l) 

write(l, ' (5x, a, f9. 5, a, a) c (7,1) ,",","56*0.0," 
write(l,'(5x,a,6(f9.5,a))')"&".c(1.2),".",c(2,2)’".".c(3,2),"," 

^  .c(4,2) , " , " ,c(5,2) . " , " ,c(6,2)  "  " 

write(l, ' (5x,a,f9.5,a,a) ')"&",c(7,2) ," ,"."56*0.0," 
write(l, '(5x,a.6(f9.5.a))')"&",c(l,3),",",c(2.3)i",".c(3,3) 

*  ,c(4,3),",",c(5.3).",",c(6,3)  "  " 

writed,  '  (5x,a,f9.5,a,a)')"&",c(7.3) ,  "  ,  "  ,  "56*0 . 0 .  " 

write (1, ' (5x,a,6(f9.5,a)) ' )"&",c(l,4) , " , " ,c(2 .4) i " , " ,c(3 ,4) . " , " 

*  ,c(4.4),",",c(5,4),",",c(6.4)  "  " 

writed,  '  (5x,  a,  £9.5,  a.  a) c(7. 4)  ,",","56*0  0/" 

write(l,'(6x,a)')"data  delta" 

write(l , ' (5x, a, 7(f8 . 3 ,a) ) ' )  "&  /" .delta(l , 1) , " , " , delta(2  1)  "  " 

&  .belta(3,l),",",delta(4.1),",",delta(5.1),",",delta(6  1)  "  " 

6c  delta(7 , 1)  , " , "  '  »  *  * 

write(l, ' (5x,a$) ')"&  " 

writed,'(7(f8.3,a))')  deltad ,  2)  ,  "  ,  "  ,  delta(2 , 2)  ,  "  " 

&  delta(3,2),".",delta(4,2),",",delta(5,2),",",dkta(6  2)  "  " 

&  delta(7 , 2) , " , "  ’  »  »  » 

writed, '  (5x,a$) ')  "&  " 

writed, '(7(f8.3,  a))')  deltad ,  3) ,  "  ,  "  ,  delta(2 . 3) ,  "  ,  "  , 

&  delta(3,3),",",delta(4,3),",",delta(5,3),",",ddta(6  3)  "  " 

&  delta(7,3),","  .  . 

writed,  '  (5x,a$)  '  )"&  " 

writed,  '(7(f8.3,  a))')  deltad  ,4) ,  "  ,  "  ,  delta(2 .4) ,  "  ,  "  , 

&  delta(3,4),",",delta(4,4),".",delta(5,4),",",ddta(6  4)  "  " 

&  delta(7 ,4) , "/"  ’  '  •  ’ 

write(l,'(6x,a)')”data  beta" 
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write(l.'(5x,a.7(f8.3.a))')  "&  /" ,beta(l . 1) . " . " .beta(2 , 1) , " , "  , 
&  beta(3.1).",".beta(4.1).".".beta(5.1) ,".",beta(6.1).".". 

&  beta(7.1) 

write(l, ' (5x,a$) ' )"&  " 

writed ,  '  (7 (f8 . 3 ,  a)  )  ' )  beta(l ,  2)  , " , "  ,beta(2 , 2) .  "  .  "  , 

&  beta(3,2).".”,beta(4,2),",",beta(5,2),".".beta(6,2) 

&  beta(7.2)."." 

write(l , ' (5x,a$) ' ) "&  " 

writed,  '(7(f8.3,a))')  betad, 3)  ,  "  ,  "  .beta(2 , 3)  .  "  .  " . 

&  beta(3 . 3) , " , " .beta(4 , 3) . " . " ,beta(5 , 3) . " , " .beta(6 , 3) . "  , "  , 

&  beta(7,3),"." 

wrlte(l,'(5x,a$)')"&  " 

write(l, '(7(f8.3,a))')  beta(l,4) . " . " .beta(2 ,4) . " . " . 

&  beta(3,4),",",beta(4,4),",",beta(5,4),",",beta(6.4) 

&  beta(7.4) 

write(l , ' (6x, a, f6 . 3 ,a) ' ) "dsta  eta  /",eta,"/" 
write(l , '  (6x,  a,  f6 . 3  ,a)  '  )"<lata  kappa  /".kappa,"/" 
writed, '  (6x, a)' )"end" 
close(l) 

return 

end 


Appendix  A  Data  Reduction/WESTEX  Texture  Prediction  Modules 


A17 


Item  A, 7:  Contents  of  file  'bdataYDP.f' 


block  data 

Bare  soil  :  Desert  Pavement,  Yuma,  AZ 

Green  vegetation:  0 
Brown  vegetation:  0 
Forest  canopy  :  0 

common  /values/  b(7 . 9 ,4) .c(7, 9 .4) ,delta(7.4) ,beta(7 ,4) ,eta,kappa 
real  eta, kappa 

data  b  /  0 . 9164613E-02 ,  0. 1064413E-01 ,  0. 1064413E-01 

6c  1,604921  ,  0.4686063E-02,  0.4686063E-02 

6c  0 . 9164613E  -02,56*0,0, 

6c  O.OOOOOOOE+00,  O.OOOOOOOE+00,  0, OOOOOOOE+00 , 

6c  O.OOOOOOOE+00,  O.OOOOOOOE+00,  0 . OOOOOOOE+00 ' 

&  0,0000000E+00, 56*0.0, 

6c  O.OOOOOOOE+OO,  0 .  OOOOOOOE+00 ,  0 ,  OOOOOOOE+00 , 

6c  0,OOOOOOOE+O0,  0 ,  OOOOOOOE+00 ,  0 ,  OOOOOOOE+00 ' 

6c  0,0000000E+00, 56*0.0, 

6c  O.OOOOOOOE+OO,  O.OOOOOOOE+OO,  O.OOOOOOOE+OO, 

6c  O.OOOOOOOE+OO,  O.OOOOOOOE+OO,  O.OOOOOOOE+OO 
6c  O.OOOOOOOE+OO,  56*0.0/ 

c  /  18.47000  ,  23,79000  ,  23.79000 

6c  107.9000  ,  36.31000  ,  36.31000 

6c  18,47000  ,56*0.0, 

6c  O.OOOOOOOE+OO,  O.OOOOOOOE+OO,  O.OOOOOOOE+OO, 

6c  O.OOOOOOOE+OO,  O.OOOOOOOE+OO,  O.OOOOOOOE+OO 
&  O.OOOOOOOE+OO, 56*0.0, 

6c  O.OOOOOOOE+OO,  O.OOOOOOOE+OO,  O.OOOOOOOE+OO 
6c  O.OOOOOOOE+OO,  O.OOOOOOOE+OO,  0 .  OOOOOOOE+00  ’ 

6c  O.OOOOOOOE+OO,  56*0.0, 

6c  O.OOOOOOOE+OO,  O.OOOOOOOE+OO,  O.OOOOOOOE+OO, 

6c  O.OOOOOOOE+OO,  O.OOOOOOOE+OO,  O.OOOOOOOE+OO 
&  O.OOOOOOOE+OO, 56*0.0/ 

data  delta  /-0.335,  0.610,  5.128,-13.22.  3.940.-4.569,-0.335 

0.000,  0.000,  0.000,  0.000,  0.000,  o.oooi  o.ooo[ 

0.000,  0.000,  0.000,  0.000,  0.000,  o.ooo]  o.ooo! 

0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000/ 

beta  /-0.127,  1.470,-1.487,-7.348,-4.205,  4.139,-0.127, 

0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000, 

0.000,  0.000,  0.000,  0.000,  0.000,  o.ooo[  o.oooi 

0.000,  0.000,  0.000,  0.000,  0.000,  0.000*  0  000/ 

data  eta  /  1.000/ 
data  kappa  /-I. 000/ 
end 
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Item  A, 8:  Contents  of  file  'bdataGTT.f' 


* 

* 

* 

* 

* 


block  data 

Bare  Soil  :  Coefficients  from  Yuma,  AZ  desert  pavement  modified 

to  simulate  test  track  at  Grayling,  MI 
Green  Vegetation:  0 
Brovm  Vegetation:  0 
Forest  Canopy  :  0 

common  /values/  b(7 , 9 ,4) ,c(7 , 9 ,4) , delta(7 , A) .beta(7 ,4) . eta.kappa 
real  eta, kappa 

data  b  /  0 . 1095171E-01 ,  0 . 1397042E-01 ,  0 . 1397042E-01 , 

&  0.4931600E+00,  0 . 1740872E-01,  0. 1740872E-01, 

&  0.1095171E-01, 56*0.0, 

&  O.OOOOOOOE+00,  O.OOOOOOOE+00,  O.OOOOOOOE+00, 

&  o.oooooooE+00,  o.oooooooE+00,  o.oooooooE+ooi 
&  O.OOOOOOOE+00, 56*0.0, 

&  O.OOOOOOOE+00,  O.OOOOOOOE+00,  O.OOOOOOOE+00, 

&  o.oooooooE+00,  o.oooooooE+00,  o.oooooooE+ooi 

&  O.OOOOOOOE+00, 56*0.0, 

&  O.OOOOOOOE+00,  O.OOOOOOOE+00,  O.OOOOOOOE+00, 

&  O.OOOOOOOE+00,  O.OOOOOOOE+00,  o.oooooooE+ooi 
6c  O.OOOOOOOE+00, 56*0.0/ 

«lata  c  /  35.54000  ,  22.70000  ,  22.70000 

6c  26.45000  ,  30.83000  ,  30.83000 

6c  35.54000  ,56*0.0, 

6c  O.OOOOOOOE+00,  0 .  OOOOOOOE+00 ,  0 .  OOOOOOOE+00 , 

6c  O.OOOOOOOE+00,  0 .  OOOOOOOE+00 ,  0 .  OOOOOOOE+00  i 

6c  O.OOOOOOOE+00, 56*0.0, 

6c  O.OOOOOOOE+00,  0 .  OOOOOOOE+00 ,  0 .  OOOOOOOE+00 , 

6c  O.OOOOOOOE+00,  0 .  OOOOOOOE+00 ,  0 .  OOOOOOOE+00 , 

&  O.OOOOOOOE+00, 56*0.0, 

6c  O.OOOOOOOE+00,  0 .  OOOOOOOE+00 ,  0.  OOOOOOOE+00 , 

6c  O.OOOOOOOE+00,  0  .  OOOOOOOE+00 ,  0 . OOOOOOOE+00 ,’ 

6c  O.OOOOOOOE+00, 56*0.0/ 

data  delta  /-0.335,  0.610,  5.128,-13.22,  3.940,-4.569,-0.335, 

6c  0.000,  0.000,  0.000,  0.000,  0.000,  o.oooi  o.ooo[ 

&  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  o.oooi 

&  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000/ 

data  beta  /-0.127,  1.470,-1.487,-7.348,-4.205,  4.139,-0.127, 

6c  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000, 

6c  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  o.ooo’ 

6c  0.000,  0.000,  0.000,  0.000,  0.000,  0.000,  0.000/ 

data  eta  /  1.000/ 
data  kappa  /-I. 000/ 
end 
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Item  A. 9:  Contents  of  file  'texgen.c' 


#include  <stdio.h> 
#include  <stdlib.h> 
#include  <string.h> 
#include  <math.h> 


typede  f  s  true  t  { 

float  atime[25], 
temp [25] , 
rad[25] , 
tempi [25] , 
radl[25], 
temp2 [25] , 
rad2[25] , 
avrad , 
minrad, 
maxrad , 
avradl , 
avrad2 ; 
int  num; 
int  day; 

)  Diurnal; 

float  ar[26]-={0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0  0  0  0  0  0  0  0  01- 

int  kflag[3]-{0,0,0) ; 

float  xlv-0.0  ,xv«-0.0,  resol; 

Diurnal  di; 
double  bandradO  ; 
float  bandradOO; 
char  str[90],  *sptr; 
int  vegindicator-0 ; 
float  sunup , sundown ; 

FILE  *in,  *texout; 
char  *clead(char  str[]) 

{ 

sptr-=str ; 

while (*sptr-”='  '  6e6£  *sptr!-0)  sptr-H-; 
return  sptr; 

) 

main (int  ac ,  char  *av[]) 

{ 

float  lemis , sabsr , sat ; 
float  thermc,  thermd; 
float  slope ; 

if(ac<6)  { 

puts("usage:  texgenOl  veg  liweg  resol  inputfile  outputfile"V 
exit(l); 
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xv-atof (av[l] ) ; 

xlv-atof (av[2] ) ; 

resol-ar[15]-atof (av[3] ) ; 

if ( (texout-fopen(av[5] , "w” ))— NULL) exit (1) ; 

if ( (in-fopen(av[4] , "rt" ))— NULL) exit (0) ; 

Se tups tandardi terns ( ) ; 

GetRads () ; 

if  (vegindicator— 1)  readvegO  ; 
else  if  (vegindicator“0)readbare()  ; 
else  if (vegindicator““2)readcan() ; 


readbareO 

( 

register  int  k,  i; 
float  Tsdv; 

float  p-  2.27E-07,  q  -  1332.5; 
float  T; 
float  TB,TB2; 


fputs ( "Bare\n" , texout) ; 

fputs  (  "Time\tSdv\tCorrL\tTemp\tBrtTmp\tSdvBT\tRad\n“ 

, texout) ; 

fputs  (*• - \t - \t . \t - \t------\t . \t - \n” 

, texout) ; 

for  (k«=0  ;k<di  ,  nujn;k++)  { 

T-di.teinp[k)+273.15; 

TB=  q  /  (logfd  +  (l/ar[3])  *  (expf(q/T)  -  1)));  /*JFMV 

ar [0]-di . rad[k] ; 

ar[l]-0.0; 

ar[2]-ar[l]; 

ar [ 16 ] “di . avrad ; 

ar[17]-0.0; 

ar[18]-ar[16] ; 

ar[19]-di.atime[k] ; 

ar ( 22 ) -ar [21] -ar [ 20 ] -di . maxrad ; 

ar [ 25 ] -ar [ 24 ] -ar [ 23 ] -di . minrad ; 
sdaclor_(&xv,&xlv,kflag, ar) ; 

Tsdv-p*TB*TB*(  exp(q/T)  -  1.0  )*ar[0] ; 

fprintf( texout , " %4 . 2f\t%7 . 2f\t%7 . 2f\t%5 . 2f\t%5 . 2f\t%7 , 2f\t%5 . 2f\n" 
,  di.atimefk] ,ar [0] ,ar[l]  , di . temp [k] ,TB-273 . 15 
,Tsdv, di . rad[k] ) ; 

) 

) 

readvegO 
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register  int  k,  i; 
float  Tsdv; 

float  p-  2.27E-07,  q  —  1332,5,  totrad; 
float  T; 
float  TB; 


SetupMVegItemsO  ; 
fputs  (  "MVeg\n"  ,  texout)  ; 

fputs ( "Time\tSdv\tCorrL\tTemp\tBrtTmp\tSdvBT\tRad\n" 

, texout) ; 

fputs("----\t - \t . \t - \t . \t . \t - \n" 

, texout) ; 

for(k-0  ;k<di  .nuin;k++)  { 

T-di.temp[k]+273.15; 

TB-  q  /  (logf(l  +  (l/ar[3])  *  (expf(q/T)  -  1)));  /*JFMV 

ar[0]-di.radl[k] ; 

ar[l]-di.rad2[k] ; 

ar[2]«ar[l] ; 

ar  [  16  ]  -=di ,  avradl ; 

ar [ 17 ] -di . avrad2 ; 

ar[18]-ar[17] ; 

ar [19]“di .atimefk] ; 

ar [ 22 ] ^ar [ 21 ] ^ar [ 20 ] -di . maxrad ; 

ar [ 2  5 ] -ar [ 24 ] -ar [23] -di . minrad ; 

sdaclor_(&xv,&xlv,kflag, ar) ; 


Tsdv=p*TB*TB*(  exp(q/T)  -  1.0  )*ar[0] ; 

fprintf( texout . "%4 . 2f\t%7 . 2f\t%7 . 2f\t%5 . 2f\t%5 . 2f\t%7 . 2f\t%5 . 2f\n" 
,  di.atime[k;]  ,ar[0]  ,ar[l]  , di . temp[k]  ,TB-273 . 15 
,Tsdv, di . rad[k] ) ; 

to trad-ar [ 9 ] *xv*di . radl [ k ] +ar [ 3 ] * ( 1 - xv) ★di . r ad2 [ k ]  ; 


readcanO 

{ 

register  int  k,  i; 
float  Tsdv; 

float  p-  2.27E-07,  q  —  1332.5,  totrad; 
float  T; 
float  TB; 


S etupCanI terns () ; 

fputs ( "Canopy\n" , texout) ; 
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fputs ( "Time\tSdv\tCorrL\tTemp\tBrtTmp\tSdvBT\tRad\n" 

, texout) ; 

fputsc - \t----\t . \t - \t . \t . \t - \n" 

, texout) ; 

for(k-0;k<di  .numjk-H-)  { 

T“di , temp[k]+273 . 15 ; 

TB-  q  /  (logf(l  +  (l/ar[3])  *  (expf(q/T)  -  1)));  /*JFM*/ 

ar[0]-di. rad[k] ; 

ar [9]-resol; 

ar[10]-di.atime[k] ; 

ar[7]-0.5; 

ar  [ 11 ] »di . avrad ; 

ar[12J-di.maxrad; 

ar [ 13 ] -di . minrad ; 


fsdaclor_(kflag,ar) ; 

Tsdv«=p*T*T*(  exp(q/T)  -  1.0  )*ar[0]  ; 

fprintf( texout . “%4 . 2f\t%7 . 2f\t%7 . 2f\t%5 . 2f\t%5 . 2f\t%7 . 2f\t%5 . 2f\n" 
,  di.atimefk] ,ar[0] ,ar[l)  ,di . temp[k] .TB-273 . 15 
,Tsdv, di. rad[k] ) ; 


} 

} 


int  GetRadsO 

{ 

float  atime , temp , rad , tempi , radl , temp 2 , rad2 ; 
register  i; 

rewind (in) ; 

di . avrad“di . avradl-di . avrad2“0 ; 
di . maxrad“0 , 0 ; 
di . minrad=200 . 0 ; 

if(vegindicator)  while (! feof( in) )  { 
if(fgets(str,80,in)— NULL)  break; 

if  (  IstrncmpC’ . REF"  ,  clead(str)  ,8)) 

/*  for  vegetation  */ 

for(i“0; i<24; i++)  { 

if(fgets(str,  80,  in)— NULL) break; 
sscanf (str , "%f  %*f  %*f  %*f  %f  %f  %f " , 
&atime ,  &temp ,  &templ , &temp2) ; 
di . atime [ i ] -atime ; 

di . temp [ i ] -temp ; 

/* 

di.rad[i]-bandrad(temp+273,15,  8.0,  12.0); 


Appendix  A  Data  Reduction/WESTEX  Texture  Prediction  Modules 


A23 


di,rad[i]-bandradO(temp+273. 15,ar[3] ) ; 

if  (di.maxrad  <  di,rad[i])  di.maxrad  -  di.rad[i] 
if  (di.minrad  >  di.rad[i))  di.minrad  -  di.rad[i] 

di . avrad  +-di . rad [ i ] ; 

di . tempi [ i ] -tempi ; 

di.radl[i]-bandrad(templ+273.15,  8.0,  12.0); 

di . radl [ i ] -bandradO ( temp 1+273 . 15 , ar [ 3 ] ) ; 
di . avradl  +-di . radl [ i ] ; 

di . temp2 [ i ] -temp2 ; 

di.rad2[i]-bandrad(temp2+273.15,  8.0,  12.0); 

di . rad2 [ i ] -bandradO ( temp2+273 . 15 , ar [ 3 ] ) ; 
di . avrad2  +-di . rad2 [ i ] ; 


if( !vegindicator)while( !feof (in))  { 
if(fgets(str,80,in)— NULL)  break; 
if ( ! strncmpC "DEG  C” , clead(str) , 5) ) 

/*  for  vegetation  */ 

for(i-0;i<24;i++)  { 

if(fgets(str,  80,  in)— NULL)break; 
s scanf ( s tr , " %  f  %  f " , &at ime , &temp )  ; 
di . atime [ i ] -atime ; 

di. temp[i]-temp; 

/* 

di.rad[i]-bandrad(temp+273.15,  8.0  12  0)* 

V 

di .  rad [ i ] -bandradO ( temp+273 . 15 , ar [ 3 ] ) ; 
if  (di.maxrad  <  di.rad[i])  di.maxrad  -  di.rad[i] 
if  (di.minrad  >  di.rad[i))  di.minrad  -  di.radfi] 
di . avrad  +-di . rad [ i ] ; 

di . tempi [ i]-0 ; 
di.radl[i]-0; 

di. temp2[i]-0; 
di.rad2[i]-0; 
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} 


di.num-i; 

if  (  di .  atime  [  i  - 1  ]  —0 . 0 )  di .  atime  [  i  - 1  ]  +-24 . 0 ; 
if  (di.num— O)return(O)  ; 
di.avrad  /-di.num; 
di . avradl  /-di . num ; 
di . avrad2  /-di . num ; 

FindSunupDownO  ; 
return(O) ; 

) 

int  FindSunupDownO 

{ 

float  atime [30],  asolar; 
register  int  i-0,  flag-0; 

rewind (in) ; 
while( ! feof (in) )  { 

if  (fgets(str ,  80 ,  in)— NULL)  break; 

if ( ! strncmp( "TIME  AIR",  clead(str),  9))break; 

) 

if (fgets(str ,  80,  in)— NULL) return  0; 

while( ! feof (in) )  { 

if (fgets (str ,  80,  in)— NULL) break; 

sscanf(str,  "%f  %*f  %*f  %*f  %*f  %f",  &atime[i],  &asolar) ; 
if(!flag  &&  asolar! -0.0)  { 

sunup— atime [ i ] ;  /*  sunu  */ 
flag-1; 

) 

else  if (flag  &&  asolar— 0.0)  { 

sundown-atime [ i-1] ;  /+  sundown  */ 
break; 

) 

i++; 

if (i>2 9) return  0; 

) 


int  Se tups tandardi terns ( ) 

{ 

float  heatcond ,  thermaldif f ,  emiss  ,  absorp ,  satur ; 
float  slope; 
register  i; 

rewind (in) ; 
while (! feof (in) )  { 
if  (  fgets  (str  ,80,  in)— NULL)  break; 
if(!strncmp( "VEGETATION",  clead(str) ,  10)) 
vegindicator— 1 ; 

if (!strncmp( "CANOPY",  clead(str) ,  6)) 
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vegindicator— 2 ; 

if(!strncmp("EMISS" ,clead(str) ,5))  ( 
fgets(str, 80, in) ; 

sscanf(str, "%f  %f  %f " ,&emiss , &absorp ,&satur) ; 

) 


if  ( !  stmcmp( "NO .  C"  ,  clead(str) ,  8)  )  { 

fgets(sti-, 80, in)  ; 

sscanf(str,"%*d  %*f  %*f  %f  %f"  ,&thennaldiff ,&heatcond) ; 


if ( !stmcmp(”SFC  SLOPE",  clead(str),  9)) 

( 

fgets(str,  80,  in); 
fgets(str,  80,  in); 
sscanf(str,  "%f",  Aslope); 

) 

) 


ar[3]-emiss; 
ar [4]-absorp; 
ar [ 5 ] -heatcond ; 
ar [ 6 ] -thermaldif f ; 
ar [ 7 ]“satur ; 
ar [ 8 ] -slope ; 

return(O) ; 

) 

int  SetupMVegItemsO 

{ 

float  folcov, emiss , absorp ; 

rewind(in) ; 
while (! feof (in) )  { 
if(fgets(str,80,in)— NULL)  break; 
i f ( I s  t rncmp ( " COVERAGE " , c lead ( s  tr ) , 8 ) )  { 
fgets(str , 80, in) ; 
fgets(str , 80 , in) ; 

sscanf (str ,  "%f  %'*’f  %f  %f",6cfolcov,  &emiss,  &absorp) ; 
break; 


) 


xv-folcov; 
ar [ 12 ] -ar [ 9 ] -emiss ; 
ar [ 13 ] — ar [ 10 ] —absorp ; 
return(O) ; 
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int  SetupCanItemsO 

{ 

float  folcov, emiss , absorp ; 
float  lai,  stoma; 

rewind  (in) ; 
while (! feof( in) )  { 
if  (fgets(str ,  80 ,  in)“NULL)  break; 
if  (  !  stmcmp( "COVERAGE"  ,  clead(str)  ,  8) )  { 
fgets(str , 80, in)  ; 
fgets(str , 80 , in) ; 

sscanf (str ,  "%f  %*f  %f  %f",6cfolcov,  &emiss,  &absorp)  ; 
break ; 

) 

) 


lai-7 . 0*folcov;  /*  this  is  strictly  made  up  */ 
stoma-. 333;  /*  so  is  this  */ 


xv-folcov; 
ar [ 2 ] — ar [ 1 ] -emiss ; 
ar [5]— absorp; 
ar [8]-stoma; 
ar [6]-lai; 

return(O) ; 

) 

float  bandradO(tmp, epsilon) 
float  tmp ; 
float  epsilon; 

{ 

float  rad; 

rad  -  3305. 8/(  expf (1332 . 5/tmp)  -  1.0  )*epsilon; 
re turn (rad) ; 

) 

/* 

double  bandrad  ( tmp ,  laml ,  lam2 ) 
double  tmp ,  laml ,  lam2 ; 

{ 

float  band_frac; 

float  next ,v,suml , sum2 ,m,mv; 

V— 14388 . 0/(laml*tmp) ; 

m-1.0; 

suml— 0 ; 

mv-m*v ; 

next—  exp(  -mv) /pow(m, 4 . 0)*(  (  (inv+3)*mv+6)*mv+6)  ; 
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do  { 

suml-h-  next; 
m++; 
mv—m^v ; 

next-  exp  (  -mv)  /pow  (m ,  4 . 0  )  *  ( (  (iiiv+3  )  *inv+6  )  *mv+6 ) ; 

)  while (  nK-lO) ; 
suml-15 . 0/97 . 40909*sLunl ; 

v-14388 .0/(lain2*tmp)  ; 

m-1.0; 

sum2-0 ; 

mv-m*v ; 

next-  exp (  -inv)/pow(m ,  4 . 0) * (  (  (mv+3 ) *iirv+6 )*iiiv+6 )  • 
do  { 

suni24—  next; 

m-H-; 

mv— m'*v; 

next-  exp ( -mv) /pow (m , 4 . 0 ) * ( ( (mv+3 ) *mv+6 ) *mv+6 ) ; 

}  while (  m<-10) ; 
sum2-15 . 0/97 . 40909*siim2 ; 

(sum2-suml)*5. 67e-8*pow(tnip,4.0))  /  3.1415926; 
return(band_frac) ; 

} 

V 


28 


Appendix  A  Data  ReductionWESTEX  Texture  Prediction  Modules 


Item  A. 10:  Contents  of  file  'sdaclor.f' 


subroutine  sdaclor(xv,xlv,kflag, a) 

****************************************************T*rit****'iHt******* 

*  This  subroutine  is  based  on  equations  26  and  27  on  page  13  of  * 

*  the  report  "Physics  Based  IR  Terrain  Radiance  Texture  Model"  * 

*  by  Richard  A.  Weiss  and  Bruce  M.  Sabol  of  the  US  Army  * 

*  Waterways  Experiment  Station,  Vicksburg,  MS  39180  for  * 

*  James  A.  Smith,  Laboratory  for  Terrestrial  Physics,  * 

*  NASA  Goddard  Space  Flight  Center,  Greenbelt,  MD  20771  * 


real  xv  !  fraction  of  terrain  covered  by  vegetation 

real  xlv  !  fraction  of  vegetation  cover  that  is  green 

real  a (22)  !  input/output  vector 
integer  kflag(3)  !  Noon/aftemoon  flag 
real  x(4) 

*********************************************************************** 


*  on  input,  a(  1)  -  radiance  of  bare  soil  * 

*  a(  2)  “  radiance  of  live  vegetation  * 

*  a(  3)  —  radiance  of  dead  vegetation  * 

*  a(  4)  “  longwave  emissivity  of  soil  * 

*  a(  5)  “  shortwave  absorptivity  of  soil  * 

*  a(  6)  -  heat  conductivity  of  soil  * 

*  a(  7)  “  thermal  diffusivity  of  soil  ★ 

*  a(  8)  -  soil  moisture  * 

*  a(  9)  -  average  terrain  slope  for  soil  * 

*  a(10)  “  longwave  emissivity  of  green  vegetation  * 

*  a(ll)  "  shortwave  absorptivity  of  green  vegetation  * 

*  ^(12)  -  average  terrain  slope  for  green  vegetation  * 

*  a (13)  -  longwave  emissivity  of  brown  vegetation  * 

*  a(14)  -«  shortwave  absorptivity  of  brown  vegetation  * 

*  a(15)  «  average  terrain  slope  for  brown  vegetation  * 

*  a(16)  -  resolution  * 

*  a(17)  -  diurnal  average  radiance  of  bare  soil  * 

*  a (18)  -  diurnal  average  radiance  of  live  vegetation  * 

*  a(19)  «=  diurnal  average  radiance  of  dead  vegetation  * 

*  a(20)  ■=  time  of  day  * 

*  a(21)  “  diurnal  maximum  radiance  of  bare  soil  * 

*  a(22)  “  diurnal  maximum  radiance  of  live  vegetation  * 

*  a(23)  -  diurnal  maximum  radiance  of  dead  vegetation  * 

*  a(24)  -  diurnal  minimum  radiance  of  bare  soil  * 

*  a (25)  -  diurnal  minimum  radiance  of  live  vegetation  * 

a(26)  -  diurnal  minimum  radiance  of  dead  vegetation  * 

*  on  return,  a(  1)  -  standard  deviation  of  radiance  * 

*  a(  2)  —  correlation  length  of  radiance  * 

*********************************************************************** 

integer  band(3)  !  band  selector 
real  u(3) ,e(3) 

real  msd(3)  !  standard  deviation  multipliers 
!  msd(l)  for  soil 

!  msd(2)  for  green  vegetation 

!  msd(3)  for  brown  vegetation 


real  msd(3) 
!  msd(l) 

!  msd(2) 

!  msd(3) 
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real  mcl(3)  !  correlation  length  multipliers 
!  mcl(l)  for  soil 

!  mcl(2)  for  green  vegetation 

!  mcl(3)  for  brown  vegetation 

real  sdvs  !  standard  deviation  of  radiance  for  'soil/vegetation 
real  civs  !  correlation  length/mixed  soil  &  vegetation 

. '  ''^^^i^nce  in  radiance  for  vegetation 

*  The  values  of  b,c, delta, beta, eta  and  kappa  are  initialized  in  * 

a  block  data  subroutine  under  the  block  "values".  In  this  way  * 

*  different  data  sets  can  be  accomodated  by  changing  out  the  * 

*  block  data  module.  ^ 

******************************^*^****^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ 
real  b(7 , 9 ,4) ,c(7 , 9 ,4)  !  vector  constants 


!  ?(band,*,l) 
!  ? (band,*, 2) 
I  ? (band,*, 3) 
!  ?(band,*,4) 


soil  vector  constants 
green  vegetation  vector  constants 
brown  vegetation  vector  constants 
forest  canopy  vector  constants 


real  delta (7, 4)  !  standard  deviation  exponents 


!  delta(band, 1) 
!  delta(band, 2) 
!  delta  (band,  3) 
!  delta(band,4) 
real  beta(7,4) 

!  beta(band,l) 

!  beta(band,2) 

!  beta (band, 3) 

!  beta(band,4) 
real  eta, kappa 


for  soil 

for  green  vegetation 
for  brown  vegetation 
for  forest  canopy 
!  correlation  length  exponents 
for  soil 

for  green  vegetation 
for  brown  vegetation 
for  forest  canopy 
!  resolution  exponents 


conunon  /values/  b(7 . 9 .4)  .c(7 , 9 ,4)  ,delta(7 .4)  ,beta(7 .4)  .eta.kappa 

**********************************^*^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ 
*  Normalize  the  radiance  values*  ^ 

*******************************^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ 
do  i=l,3 

if  (  a(i+16)  .eq.  0.0  )  then 
a(i)  -  0.0 
else 

a(i)  -  a(i)/a(i+16) 
end  if 
end  do 


*  Determine  which  band  to  retrieve  values  from:  * 

do  i-1,3 

x(l)  —  a(i)  !  radiance 

x(2)  —  a(i+16)  !  diurnal  average  radiance 

x(3)  “  a(i+20)  !  diurnal  maximum  radiance 

x(4)  -=  a(i+23)  !  diurnal  minimum  radiance 

call  adj (x,b, delta, i) 
call  adj (x,c, beta, i) 
call  ge t__band(i,kf lag, X, band) 
end  do 
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St********************************************************************** 

*  Compute  standard  deviation  and  correlation  length  multipliers  * 

*  for  bare  soil.  * 

*********************************************************************** 

msd(l)  -  0.0 
mcl(l)  -  0.0 
do  i  -  4,9 

msd(l)  -  msd(l)  +  b(band(l) , i-2 , l)*a(i) 
mcl(l)  -  mcl(l)  +  c(band(l) ,i-2,l)*a(i) 
end  do 

msd(l)  -  msd(l)  +  b(band(l) , 1 , 1) 
mcl(l)  -  mcl(l)  +  c(band(l),l,l) 

*  Compute  standard  deviation  and  correlation  length  multipliers  * 

*  for  green  vegetation.  * 

*********************************************************************** 

msd(2)  -=  0.0 
mcl(2)  -  0.0 
do  i  -=  10,12 

msd(2)  msd(2)  +  b(band(2)  ,  i-8 , 2)*a(i) 
mcl(2)  «  mcl(2)  +  c (band(2) , i-8 , 2)*a(i) 
end  do 

msd(2)  -  msd(2)  +  b(band(2) , 1 , 2) 
mcl(2)  -  mcl(2)  +  c(band(2) , 1 , 2) 

****************-*-***********TSr************************-)»:*****‘)V*****:*:***** 

*  Compute  standard  deviation  and  correlation  length  multipliers  * 

*  for  brown  vegetation.  * 

msd(3)  —  0.0 
mcl(3)  =  0.0 
do  i  -=  13,15 

msd(3)  -=  msd(3)  +  b(band(3)  ,  i-ll,3)*a(i) 
mcl(3)  -  mcl(3)  +  c(band(3) , i-11 , 3)*a(i) 
end  do 

msd(3)  -  msd(3)  +  b(band(3) , 1 , 3) 
mcl(3)  -  mcl(3)  +  c(band(3) , 1 , 3) 


Equation  26,  page  13  of  report  cited  above: 


*  on  input,  a(  1)  —  radiance  of  bare  soil  * 

*  a(  2)  “  radiance  of  live  vegetation  * 

*  a(  3)  -  radiance  of  dead  vegetation  * 

*  a(16)  “  resolution  * 

***************************************************Tt***-*-*****-*-*******-3^ 

do  i-=l ,  3 

if  (  a(i)  .le.  0.0  )  then 
u(i)  -  0.0 
else 

u(i)  -=  a(i)**delta(band(i)  ,  i) 
end  if 
end  do 
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varv  -  (1.0  -  xlv)*(msd(3)*u(3))**2  +  xlv*(msd(2)*u(2) )**2 

sdvs  -  a(16)**eta*sqrt((1.0-xv)*(msd(l)*u(l))**2+xv*varv)  - 

***************************************^*^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ 
*  Equation  27,  page  13  of  report  cited  above-  * 

**************************************^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ 


do  i-1,3 

if  (  a(i)  .le.  0.0  )  then 

V  -  0.0 

else 

V  -  a(16)**kappa*iiicl(i)*a(i)**beta(band(i)  i) 

end  if 

if  (  V  .eq.  0.0  )  then 
e(i)  -  0.0 
else 

e(i)  -  exp(-1.0/v) 
end  if 
end  do 

r  -  (1.0-xv)*e(l)  +  xv*(  (1 . 0-xlv)*e(3)  +  xlv*e(2)  ) 
civs  -  -1.0/alog(r) 

a(l)  —  sdvs 
a(2)  -=  civs 

return 

end  !  sdaclor 


subroutine  get_band( i , kf lag , xx . band) 

real  xx(4) 
integer  band(3) 
integer  kflag(3) 

data  s  /0.5/ 
data  t  /0.84/ 

X  -=  xx(l) 
rn  -  xx(4)/xx(2) 
rx  -  xx(3)/xx(2) 
ra  =  (1.0  -  s)*rn  +  s 
rb  -  1.0  -  t  +  t*rx 


if  (  kflag(i)  .eq.  0  .and.  x  .It.  ra  )  then 
band(i)  -  1 

°  -1®-  )  then 

band(i)  -=  4 

if  (  rx  -  X  .le.  0.0001  )  then 
kflag(i)  -  1 
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end  if 

else  if  (  kflag(i)  .eq.  1  .and.  1.0  ,le.  x  )  then 
band(i)  -=  5 

else  if  (  kflag(i)  .eq.  1  .and.  ra  .le.  x  .and.  x  .It.  1.0  )  then 
band(i)  ■=  6 
else 

band(i)  -  7 
end  if 


return 

end 

************************************************************************ 
subroutine  adj(x,a,b,i) 
real  x(4) , a(7 , 9 ,4) ,b(7 , 9 ,4) 
data  s  /0,50/ 
data  t  /0.84/ 

rmax  -  x(3)/x(2) 

rmin  -  x(4)/x(2) 

ra  =  (1.0  -  s)*rmin  +  s 

rb  -  1.0  -  t  +  t*rmax 

al  «  a(l,l,i) 

a3  -  a(3,l,i) 

a5  -  a(5,l.i) 

bl  =  b(l,l,i) 

b3  "  b(3,l,i) 

b5  -  b(5,l,i) 


a(2,l.i) 

a(6,l.i) 

b(2,l.i) 

b(4,l,i) 

b(6,l.i) 

a(4.1,i) 


a3 

a5 

bl  +  (alog(al) -alog(a3) )/alog(ra) 
alog(  a3*rb**b3/(a5*rmax**b5)  )/alog(rb/rmax) 
bl  +  (alog(al)  -  alog(a5) )/alog(ra) 
a5*rmax**(b5  -  b(4,l,i)) 


return 

end 
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Item  A. 11:  Contents  of  file  ' fsdaclor . f ' 


subroutine  fsdaclor (kflag, a) 

*  This  subroutine  is  based  on  equations  50  and  51  on  page  18  of  * 

*  the  report  "Physics  Based  IR  Terrain  Radiance  Texture  Model"  * 

*  by  Richard  A.  Weiss  and  Bruce  M.  Sabol  of  the  US  Army  * 

*  Waterways  Experiment  Station,  Vicksburg,  MS  39180  for  ★ 

*  James  A.  Smith,  Laboratory  for  Terrestrial  Physics,  * 

*  NASA  Goddard  Space  Flight  Center,  Greenbelt  MD  20771  * 

real  a (14)  !  input/output  vector 

*  on  input,  a(  1)  -  computed  radiance  of  forest  canopy  * 

a(  2)  *  longwave  emissivity  of  leaves  * 

a(  3)  -  longwave  absorptivity  of  leaves  * 

*  a(  4)  -  longwave  emissivity  of  soil  * 

a(  5)  -=  longwave  absorptivity  of  soil  * 

a(  6)  -  shortwave  absorptivity  of  leaves  * 

*  a(  7)  ■«  leaf  area  index  * 

*  a(  8)  -  leaf  slope  distribution  * 

*  a(  9)  -  stomatal  resistance  to  water  vapor  diffusion  * 

a(10)  »=  resolution  * 

*  a(ll)  -=  time  of  day  * 

a (12)  *=  average  diurnal  radiance  of  canopy 

a (13)  =  maximum  diurnal  radiance  of  canopy  * 

a (14)  -  minimum  diurnal  radiance  of  canopy  * 

*  on  return,  a(  1)  «  standard  deviation  of  radiance  * 

*  a(  2)  -  correlation  length  of  radiance  * 

real  x(4) 


real  b(7 , 9 ,4) ,c(7 , 9 ,4) 
!  ?(band,*,l) 

!  ?(band,*,2) 

!  ? (band,*, 3) 

!  ? (band,*, 4) 
real  delta(7,4)  !  stan 
!  delta(band, 1)  for 

!  delta(band, 2)  for 

!  delta (band , 3 )  for 

!  delta(band,4)  for 

real  beta(7,4)  !  corr 

!  beta(band,l)  for 

!  beta (band, 2)  for 

!  beta (band, 3)  for 

!  beta (band, 4)  for 


9,4)  !  vector  constants 

soil  vector  constants 
green  vegetation  vector  constants 
brown  vegetation  vector  constants 
forest  canopy  vector  constants 
standard  deviation  exponents 
for  soil 

for  green  vegetation 
for  brown  vegetation 
for  forest  canopy 
correlation  length  exponents 
for  soil 


!  beta(band,2)  for  green  vegetation 

!  beta(band,3)  for  brown  vegetation 

!  beta (band, 4)  for  forest  canopy 

real  eta, kappa  !  resolution  exponents 

common  /values/  b(7 , 9 . 4) , c(7 , 9 ,4) , delta(7 , 4) ,beta(7 .4)  eta  kanna 

'V'o  a  T  I  ^ j 1 


real  msd 
real  mcl 
real  sd 


standard  deviation  multiplier 
correlation  length  multiplier 
standard  deviation  of  radiance 
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real  cl  !  correlation  length  of  radiance 

integer  band  !  time  of  day  band 


*  Normalize  radiance 

if  (  a(12)  .le.  0.0  )  then 
a(l)  -  0.0 
else 

a(l)  -  a(l)/a(12) 
end  if 


x(l)  -  a(l) 
x(2)  -  a(12) 
x(3)  -  a(13) 
x(4)  -  a(14) 

call  adj (x,b,delta,4) 

call  adj (x, c ,beta ,4) 

call  get_band(l ,kf lag, X, band) 

*********************************************************************** 
*  Compute  standard  deviation  and  correlation  length  multipliers  * 


msd  —  0.0 
mcl  “  0.0 
do  i  -  2,9 

msd  —  msd  +  b(band, i-1 ,4)*a(i) 
mcl  —  mcl  +  c(band, i-l,4)*a(i) 
end  do 

msd  —  msd  +  b (band, 1,4) 
mcl  —  mcl  +  c (band, 1,4) 


*  Equation  54,  55,  pages  19,20  of  report  cited  above:  * 

*  on  input,  a(  1)  —  radiance  of  bare  soil  * 

*  a(10)  -  resolution  * 

★*************************************************************T*:******** 

if  (  a(l)  .le,  0.0  )  then 
sd  -=  0.0 
cl  -  0.0 
else 

sd  -*  a(10)**eta*a(l)**delta(band,4)*msd 
cl  —  a(10)**kappa*a(l)**beta(band,4)*mcl 
end  if 

a(l)  -  sd 
a(2)  -  cl 

return 

end  !  fsdaclor 
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Appendix  B 
Data  Files 


Appendix  B  Data  Files 


Item  B.l:  Contents  of  file  'ss040893.dat'  (Yuma,  AZ) 


Resolution  Range:  0.01853  0.01668 

FEATURE  IMAGE  STD.  CORRELATON 

DATETIME  DESCRIPTION  MEAN  DEV.  LENGTH 

W/(M'2*SR)  W/(M^2*SR)  (METERS) 

08APR93:00:34  GRASS -BARESOIL  28.575  0.5002  0.3100 

08APR93 : 01 : 04  GRASS-BARESOIL  29.182  0.4761  0^3000 

08APR93:01:13  GRASS-BARESOIL  29.255  0.4686  0.3600 

08APR93.02:39  GRASS-BARESOIL  28.313  0.5103  0.3000 

08APR93:03:14  GRASS-BARESOIL  27.703  0.4443  0^2600 

08APR93:03:20  GRASS-BARESOIL  27.740  0.4452  0.2200 

08APR93:04:16  GRASS-BARESOIL  26.441  0.3373  0  1800 

08APR93: 04:40  GRASS-BARESOIL  26.827  0.3571  0^2100 

08APR93:05:05  GRASS-BARESOIL  26.827  0.4202  o!2800 

08APR93:05:25  GRASS-BARESOIL  27.971  0.5620  o!2900 

08APR93:06:13  GRASS-BARESOIL  27.896  0.5018  0^2800 

08APR93:06:32  GRASS-BARESOIL  28.009  0.4683  0^2800 

08APR93:07:12  GRASS-BARESOIL  31.092  0.6329  0.3100 

08APR93:07:28  GRASS-BARESOIL  32.943  0.5968  0  3400 

08APR93:08:07  GRASS-BARESOIL  35.172  0.4137  0  4500 

08APR93:08:28  GRASS-BARESOIL  36.723  0.3555  0  4200 

08APR93:09:16  GRASS-BARESOIL  39.793  0.4108  0^3600 

08APR93:09:32  GRASS-BARESOIL  40.677  0.5411  0  4000 

08APR93:10:10  GRASS-BARESOIL  42.947  0.7266  0  3000 

08APR93:10:29  GRASS-BARESOIL  43.820  0.8384  0  2700 

08APR93:11:06  GRASS-BARESOIL  45.024  0.8719  0  2900 

08APR93:11:31  GRASS-BARESOIL  45.951  1.0743  0  2400 

08APR93:12:08  GRASS-BARESOIL  47.594  1.1358  0  1600 

08APR93:12:32  GRASS-BARESOIL  47.492  1.3077  0  2000 

08APR93:13:04  GRASS-BARESOIL  48.305  1.3966  0  2500 

08APR93:13:19  GRASS-BARESOIL  48.170  1.2517  o!l600 

08APR93:14:14  GRASS-BARESOIL  50.432  1.3791  0  2200 

08APR93:14:27  GRASS-BARESOIL  49.796  1.4170  0  1400 

08APR93:15:09  GRASS-BARESOIL  49.927  1.1929  o!l500 

08APR93:15:34  GRASS-BARESOIL  48.543  1.0595  o!llOO 

08APR93:16:03  GRASS-BARESOIL  48.067  1.0793  0  1500 

08APR93:16:18  GRASS-BARESOIL  47.320  0.8920  0  1500 

08APR93:18:09  GRASS-BARESOIL  41.904  0.2828  0  2000 

08APR93:18:28  GRASS-BARESOIL  40.572  0.3489  0^2900 

08APR93:19:07  GRASS-BARESOIL  37.691  0.3938  0.6100 

08APR93:19:32  GRASS-BARESOIL  36.041  0.2935  0  3300 

08APR93:20:03  GRASS-BARESOIL  34.699  0.2445  0^2900 

08APR93:20:29  GRASS-BARESOIL  33.528  0.2622  0^2400 

08APR93:21:05  GRASS-BARESOIL  32.870  0.2572  0^3400 

08APR93:21:34  GRASS-BARESOIL  32.020  0.2899  0  3800 

08APR93:22:06  GRASS-BARESOIL  31.536  0.3586  0^2500 

08APR93:22:36  GRASS-BARESOIL  31.203  0.4163  0^2700 

08APR93:23:10  GRASS-BARESOIL  31.203  0.4482  0  2200 

08APR93:23:22  GRASS-BARESOIL  30.832  0.3793  0.2200 
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Item  B.2:  Contents  of  file  'ss042693.dat'  (Yuma,  AZ) 


Resolution  Range 
IMAGE 

DATETIME 

:  0.02363  0, 

FEATURE 
DESCRIPTION 

.02041 

IMAGE 

MEAN 

W/(M"2*SR) 

STD. 

DEV. 

W/(M"2*SR) 

CORRELATON 

LENGTH 

(METERS) 

26APR93:00:06 

SLOPE -DESPAV 

31.908 

0.3989 

0.4400 

26APR93:00:20 

SLOPE -DESPAV 

32.168 

0.4600 

0.4700 

26APR93:01;11 

SLOPE-DESPAV 

32.503 

0.3898 

0.4000 

26APR93:01:36 

SLOPE-DESPAV 

32.647 

0.5580 

0.3900 

26APR93:02:03 

SLOPE-DESPAV 

32.539 

0.5141 

0.3800 

26APR93:02:23 

SLOPE-DESPAV 

32.575 

0.6027 

0.3900 

26APR93:03:04 

SLOPE-DESPAV 

31.760 

0.4074 

0.4000 

26APR93:03:36 

SLOPE-DESPAV 

31.057 

0.4023 

0.4500 

26APR93:04:02 

SLOPE-DESPAV 

31.057 

0.4569 

0.4800 

26APR93:04:33 

SLOPE-DESPAV 

30.383 

0.3762 

0.4600 

26APR93:05:03 

SLOPE-DESPAV 

30.235 

0.4456 

0.4600 

26APR93:05:13 

SLOPE-DESPAV 

30.195 

0.4531 

0.4500 

26APR93:06:03 

SLOPE-DESPAV 

30.195 

0.4639 

0.3800 

26APR93:06:14 

SLOPE-DESPAV 

30.757 

0.5331 

0.3400 

26APR93:07:02 

SLOPE-DESPAV 

34.774 

0.6056 

0.3000 

26APR93:07:33 

SLOPE-DESPAV 

37.837 

0.4779 

0.3000 

26APR93:08:08 

SLOPE-DESPAV 

41.625 

0.4227 

0.3600 

26APR93:08:21 

SLOPE-DESPAV 

43.991 

0.5402 

0.4800 

26APR93:09:01 

SLOPE-DESPAV 

47.085 

0.8603 

0.4800 

26APR93:09:48 

SLOPE-DESPAV 

49.626 

1.3169 

0.3700 

26APR93:10:01 

SLOPE-DESPAV 

51.033 

1.4664 

0.4300 

26APR93:10:27 

SLOPE-DESPAV 

53.664 

1.6298 

0.4000 

26APR93:11:05 

SLOPE-DESPAV 

54.951 

1.7390 

0.4300 

26APR93:11:33 

SLOPE-DESPAV 

55.774 

1.8408 

0.4200 

26APR93:12:01 

SLOPE-DESPAV 

55.774 

1.6861 

0.3100 

26APR93:12:18 

SLOPE-DESPAV 

56.491 

1.5996 

0.3600 

26APR93:13:02 

SLOPE-DESPAV 

57.509 

1.1794 

0.2500 

26APR93:13:26 

SLOPE-DESPAV 

57.178 

1.0663 

0.2800 

26APR93:14:08 

SLOPE-DESPAV 

56.981 

1.0279 

0.2800 

26APR93:14:43 

SLOPE-DESPAV 

56.596 

0.7786 

0.2900 

26APR93:15:08 

SLOPE-DESPAV 

56.198 

0.6486 

0.2000 

26APR93:15:19 

SLOPE-DESPAV 

56.393 

0.6901 

0.1900 

26APR93:16:03 

SLOPE-DESPAV 

54.689 

0.4507 

0.2800 

26APR93:16:35 

SLOPE-DESPAV 

52.340 

0.4064 

0.2400 

26APR93:17:11 

SLOPE-DESPAV 

51.971 

0.3904 

0.5400 

26APR93:17;22 

SLOPE-DESPAV 

50.366 

0.3738 

0.5100 

26APR93:18:04 

SLOPE-DESPAV 

47.047 

0.2939 

0.7500 

26APR93:18:19 

SLOPE-DESPAV 

46.089 

0.2834 

0.5800 

26APR93:19:08 

SLOPE-DESPAV 

42.220 

0.3900 

0.5600 

26APR93:19:28 

SLOPE-DESPAV 

40.642 

0.3151 

0.4700 

26APR93:20:03 

SLOPE-DESPAV 

38.942 

0.3592 

0.5000 

26APR93:20:32 

SLOPE-DESPAV 

37.911 

0.2940 

0.5000 

26APR93:21:03 

SLOPE-DESPAV 

37.194 

0.3398 

0.5400 

26APR93:21:32 

SLOPE-DESPAV 

36.690 

0.3018 

0.5400 

26APR93:22:05 

SLOPE-DESPAV 

36.579 

0.3324 

0.5200 

26APR93:22:33 

SLOPE-DESPAV 

36.403 

0.5190 

0.3900 
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B3 


26APR93 :23 :03  SLOPE-DESPAV 
26APR93 :23: 34  SLOPE-DESPAV 


35.643  0.2804 

35.318  0.4922 


0.3900 

0.3700 
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Item  B.3:  Contents  of  file  ' ss040893 .prm' 


eta  -  1.00 

kappa  —  *1.00 

08APR9  3  SLOPE - DESPAV 

Bandsize  1,3, 4, 5:  16  6  5  7 

Minimum  Radiance  /  Average  Radiance  -  0.705 

Maximum  Radiance  /  Average  Radiance  -  1.404 


band 

b 

delta 

c 

beta 

1 

0.8540429E-02 

-0.618 

27.10 

1.081 

2 

0.1022623E-01 

0.510 

25.87 

0.790 

3 

0.1022623E-01 

4.963 

25.87 

-2.255 

4 

1.924964 

-13.484 

39.40 

-4.011 

5 

0.7756619E-02 

1.430 

23.83 

-2.831 

6 

0.7756619E-02 

-1.221 

23.83 

0.275 

7 

0.8540429E-02 

-0.618 

27.10 

1.081 
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Item  B.4:  Contents  of  file  ' ss042693 .prm' 

eta  -  1.00 

kappa  -  -1.00 

26APR93  SLOPE-DESPAV 

Bands ize  1,3, 4, 5:  17  6  4  11 

Minimum  Radiance  /  Average  Radiance  -  0.709 

Maximum  Radiance  /  Average  Radiance  -  1,351 


band 

b 

delta 

c 

beta 

1 

0.9793684E-02 

-0.069 

12,88 

-1.264 

2 

0.1107910E-01 

0,716 

21.88 

2,111 

3 

0.1107910E-01 

5.294 

21.88 

-0.718 

4 

1.278622 

-12.884 

380.12 

-11.520 

5 

0.3400406E-02 

5.538 

47.47 

-5.080 

6 

0.3400406E-02 

-6.803 

47.47 

7.042 

7 

0.9793684E-02 

-0.069 

12.88 

-1.264 

B6 
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Item  B.5:  Contents  of  file  'ss.all' 


eta  -  1.00 

kappa  -  -1.00 

08APR93  SLOPE-DESPAV 

Bandsize  1,3, 4, 5:  16  6  5  7 

Minimum  Radiance  /  Average  Radiance  -  0.705 

Maximum  Radiance  /  Average  Radiance  -  1.404 


band 

b 

delta 

c 

beta 

1 

0.8540429E-02 

-0.618 

27.10 

1.081 

2 

0.1022623E-01 

0.510 

25.87 

0.790 

3 

0.1022623E-01 

4.963 

25.87 

-2.255 

4 

1.924964 

-13.484 

39.40 

-4.011 

5 

0.7756619E-02 

1.430 

23.83 

-2.831 

6 

0.7756619E-02 

-1.221 

23.83 

0.275 

7 

0.8540429E-02 

-0.618 

27.10 

1.081 

26APR93  SLOPE-DESPAV 

Bandsize  1, 3,4,5:  17 

6  4 

11 

Minimum  Radiance  /  Average  Radiance  -  0.709 
Maximum  Radiance  /  Average  Radiance  -  1.351 


band 

b 

delta 

c 

beta 

1 

0.9793684E-02 

-0.069 

12.88 

-1.264 

2 

0.1107910E-01 

0.716 

21.88 

2.111 

3 

0.1107910E-01 

5.294 

21.88 

-0.718 

4 

1.278622 

-12.884 

380.12 

-11.520 

5 

0.3400406E-02 

5.538 

47.47 

-5.080 

6 

0.3400406E-02 

-6.803 

47.47 

7.042 

7 

0.9793684E-02 

-0.069 

12.88 

-1.264 
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Item  B.7:  Contents  of  file  'model. db'  (same  as  modelYDP.db) 


eta  -  1.00 
kappa  -  -1.00 
BARE  SOIL 


band 

b 

delta 

c 

beta 

1 

0.9164613E-02 

-0.335 

18.47 

-0.127 

2 

0.1064413E-01 

0.610 

23.79 

1.470 

3 

0.1064413E-01 

5.128 

23.79 

-1.487 

4 

1.604921 

-13.217 

107.90 

-7.348 

5 

0,4686063E-02 

3.940 

36.31 

-4.205 

6 

0.4686063E-02 

-4.569 

36.31 

4.139 

7 

G.9164613E-02 

-0.335 

18.47 

-0.127 

GREEN 

VEGETATION 

band 

b 

delta 

c 

beta 

1 

O.OOOOOOOE+00 

0.000 

0.00 

0.000 

2 

O.OOOOOOOE+00 

0.000 

0.00 

0.000 

3 

O.OOOOOOOE+OO 

0.000 

0.00 

0.000 

4 

O.OOOOOOOE+00 

0.000 

0.00 

0.000 

5 

O.OOOOOOOE+OO 

0.000 

0.00 

0.000 

6 

O.OOOOOOOE+OO 

0.000 

0.00 

0.000 

7 

O.OOOOOOOE+OO 

0.000 

0.00 

0.000 

BROWN 

VEGETATION 

band 

b 

delta 

c 

beta 

1 

O.OOOOOOOE+OO 

0.000 

0.00 

0.000 

2 

O.OOOOOOOE+OO 

0.000 

0.00 

0.000 

3 

O.OOOOOOOE+OO 

0.000 

0.00 

0.000 

4 

O.OOOOOOOE+OO 

0.000 

0.00 

0.000 

5 

O.OOOOOOOE+OO 

0.000 

0.00 

0.000 

6 

O.OOOOOOOE+OO 

0.000 

0.00 

0.000 

7 

O.OOOOOOOE+OO 

0.000 

0.00 

0,000 

FOREST  CANOPY 

band 

b 

delta 

c 

beta 

1 

O.OOOOOOOE+OO 

0.000 

0.00 

0.000 

2 

O.OOOOOOOE+OO 

0.000 

0.00 

0.000 

3 

O.OOOOOOOE+OO 

0.000 

0.00 

0.000 

4 

O.OOOOOOOE+OO 

0.000 

0.00 

0.000 

5 

O.OOOOOOOE+OO 

0.000 

0.00 

0.000 

6 

O.OOOOOOOE+OO 

0.000 

0.00 

0.000 

7 

O.OOOOOOOE+OO 

0.000 

0.00 

0.000 
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